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ABSTRACT 


This  study  presents  new  analytical  findings  and  a 
new  computational  method  for  time-scale  decoupling  and  order 
reduction  of  linear  control  systems.  A transformation  of 
variables  is  considered  which  reduces  the  system  to  block- 
diagonal  form  with  slow  and  fast  modes  decoupled.  Once  the 
decoupled  form  of  the  system  is  obtained,  a small  parameter 
can  be  identified  which  provides  a measure  of  the  system's 
time  scale  separation,  and  the  methods  of  singular  perturba- 
tions can  be  applied  to  yield  reduced-order  approximations  of 
the  original  system.  A particular  reduced-order  model  so 
obtained  overcomes  certain  difficulties  encountered  with 
other  well  known  order-reduction  methodologies. 

The  type  of  control  system  considered  here  is  a very 
general,  widely  used  form.  The  ability  to  identify  a 
specific  small  parameter  in  all  such  systems  provides  a 
useful  link  to  singular  perturbations  theory  which  requires 
a small  parameter  to  be  identified.  Although  the  examples 
of  two-time-scale  systems  presented  here  are  relatively  low 
order,  due  to  recent  advances  in  computational  methods  for 
matrices  these  decoupling  techniques  can  be  directly  applied 
to  large  scale  systems  with  orders  of  100  or  more. 

ix 


X 


In  Chapter  I,  a two-time-scale  property  is  defined  in 
terms  of  system  eigenvalues,  and  a general  class  of  linear 
transformations  which  yield  time  scale  decoupling  is  also 
defined. 

Chapter  2 considers  two  types  of  decoupling  trans- 
formations, classical  modal  matrix  transformations  which 
decouple  all  modes  of  the  system,  and  a new  block-decoupling 
transformation  introduced  in  a short  1975  paper  by  Petar  V. 
Kokotovic  which  appeared  in  the  IEEE  Transactions  on  Auto- 
matic Control.  Matrices  L and  K used  in  this  transforma- 
tion satisfy  a nonsymmetric  algebraic  Riccati  equation  and  a 
Lyapunov  equation,  respectively. 

Application  of  the  block-decoupling  transformation  to 
reduced  order  modeling  is  considered.  Past  work  on  reduced 
order  modeling  based  on  modal  analysis  is  described  and  a 
new  model  based  on  a zeroth  order  outer  solution  is  intro- 
duced. The  differences  in  dynamic  characteristics  between 
new  and  old  reduced  order  models  are  illustrated  with  a 
second-order  example.  Feedback  matrix  design  using  block 
decoupling  transformations  is  described. 

In  Chapter  3,  the  nonsymmetric  algebraic  Riccati 
equation  is  studied  and  a basic  theorem  establishing  new 
conditions  for  existence  of  solutions  is  presented.  It  is 
proven  that  for  any  given  two-time-scale  system  the  L and 
K matrices  used  in  the  decoupling  transformation  are  real 


and  unique.  The  resulting  decoupling  transformation  is 
always  nonsingular. 

Chapter  4 presents  a new  method  of  computing  the  L 
and  K matrices.  It  is  shown  that  the  method  has  an 
asymptotic  rate  of  convergence  equal  to  the  system's  small 
parameter,  and  therefore  converges  rapidly  for  strongly  two- 
time-scale  systems.  The  method  is  initialized  using  certain 
of  the  system  eigenvectors  corresponding  to  either  the  system 
slow  or  fast  modes,  and  then  highly  accurate  L and  K 
matrices  can  be  obtained  by  linear  interation. 

In  Chapter  4,  the  choice  of  the  order  of  the  slow  sub- 
system and  the  ordering  of  state  variables  is  discussed.  The 
relationship  between  ordering  and  the  norms  of  the  L and  K 
matrices  is  described  using  specific  slow  and  fast  mode 
coupling  ratios. 

The  stability  and  rate  of  convergence  of  the  numerical 
decoupling  method  and  the  accuracy  of  the  new  reduced  order 
model  is  demonstrated  in  Chapter  5 using  two  examples  from 
recent  literature.  They  are  a fourth  order  model  of  F8  air- 
craft longitudinal  dynamics  and  a sixteenth  order  model  of 
the  F100  turbofan  engine. 

The  time-scale  decoupling  transformation  is  extended 
in  Chapter  6 to  time-varying  linear  systems.  Methods  of  com- 
puting a time-varying  LK  transformation  are  presented  and 
illustrated  with  a second-order  example. 
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CHAPTER  1 

INTRODUCTION 

Physical  systems  often  exhibit  behavior  that  includes  i 

\ 

both  fast  and  slow  dynamic  response,  and  it  is  common  engi- 
neering practice  to  separately  model  and  study  the  fast  and 
slow  phenomena.  Such  dynamic  separation  can  simplify  both  1 

analytical  and  computational  tasks,  and  becomes  a necessity  1 

for  the  study  of  large-scale  (i.e.,  large  order)  systems,  I 

such  as  those  encountered  in  structural  dynamics  (Foss,  1958) , i 

and  power  distribution  systems  (Van  Ness,  1977)  . Past  ap- 
proaches to  the  separation  of  fast  and  slow  dynamics  have 
either  relied  upon  engineering  judgment  to  choose  some  vari- 
ables as  slow  (and  neglect  the  exact  modal  structure  of  the  j , 

system) , or  employed  modal  variables  which  lack  physical 
meaning  and  are  computationally  awkward.  Here  we  present  a 
new  approach  to  the  analysis  of  such  systems  which  both 
clarifies  the  fast/slow  dynamic  separation  task  and  overcomes 
common  computational  difficulties. 

This  study  will  be  limited  to  dynamic  systems  which 
can  be  approximated  by  a set  of  linear  first-order  dif- 
ferential equations  of  the  form 

^ = x = Ax  + B(t)u,  0 < t < T (1.1) 

r 


i 


2 


where  x is  an  n-dimens-'  -'nal  vector  of  state  variables 
with  x(0)  often  specified,  u = u(t)  is  an  m-dimensional 
vector  of  control  variables  generally  specified  on  [0,T], 
and  A and  B(t)  are  real  matrices  of  appropriate  dimen- 
sions. Such  systems  often  include  a p-dimensional  vector 
of  output  variables  w,  which  is  a linear  combination  of 
the  state  and  control 


w = C ( t)  x + D (t)  u , (1.2) 

C and  D being  real.  The  matrices  B,  C,  and  D are 
assumed  to  be  continuous  functions  of  time  t on  [0,T] 
and  A is  assumed  to  be  constant.  The  systems  (1.1),  (1.2) 
will  be  referred  to  as  a control  system,  and  the  case  of 
control  systems  with  time  varying  A matrix  is  considered 
in  a later  chapter.  In  applications,  the  variables  x,  u 
and  w will  usually  represent  physical  quantities. 

As  proposed  by  Chow  and  Kokotovic  (1976) , the  control 
system  will  be  called  strongly  two-time-scale  or  simply  two- 
time-scale  if  the  eigenvalues  of  A,  represented  as  A (A) , 
can  be  separated  by  absolute  values  into  nonempty  sets  S 
and  F with  nn  elements  s.  and  n~  elements  f.  re- 
spectively,  such  that 


| 8 i | <<  | f j | for  all  s.  in  S and  fj  in  F, 


(1.3) 


I 


3 

and  = n.  If  this  relationship  is  only  satisfied 

for  the  weak  inequality 


s . < f . 

i1  j 


for  all  s.  in  S and  f.  in  F, 
i 3 


(1.4) 


the.,  the  control  system  will  be  referred  to  as  weakly  two- 
time-scale.  The  special  case  of  control  systems  which  do 
not  satisfy  (1.4),  i.e.,  for  which  all  eigenvalues  of  A 
have  equal  modulus,  will  only  be  treated  briefly  in  a later 
chapter.  This  algebraic  definition  of  two-time-scale  con- 
trol systems  does  not  necessarily  imply  a clearly  identi- 
fiable slow/fast  structure  in  the  dynamic  response  of  x(t) 
or  w(t)  which  depend  upon  the  specific  control  matrix  B(t) 
and  control  input  u(t). 

The  reason  for  using  the  slow/fast  terminology  be- 
comes clear  if  one  considers  the  homogeneous  problem  with 
zero  control.  Likewise,  for  problems  with  slowly  varying 
B ( t) u ( t)  and  for  time-intervals  longer  than  the  fast  mode 
time  constant  Minl'RefTJ'  one  can  exPect  t*ie  dynamic  re- 
sponse to  often  have  a boundary  layer  structure  at  the  end- 
points corresponding  to  the  fast  eigenvalues  f ^ plus  a 
slowly- varying  part  throughout  [0,T]  corresponding  to  the 
slow  eigenvalues  s^. 

For  any  two-time-scale  system  we  can  order  the 


eigenvalues  by  modulus  so  that 


4 


1 <_  i <_  n1  - 1 


f . 
3 


i lfj+1l 


1 <_  j <_  n2  - 1 


then  the  ratio 

V - Is  | / |f,|  « 1 d.5) 

nl  x 

will  define  the  system  small  parameter  which  provides  a 
measure  of  the  system's  time  scale  separation  and  identifies 
the  system  as  a singular  perturbations  problem.  Such  systems 
with  slow  and  fast  variables  are  also  described  in  the  liter- 
ature as  stiff  (Flaherty  and  O'Malley  1977,  and  Willoughby 
1974)  or  illconditioned  (Kokotovic  1975).  This  particular 
choice  of  small  parameter  was  also  proposed  by  Kelley  (1973). 

Other  additional  constraints  are  often  placed  on  the 
spectrum  of  the  A matrix  to  insure  stability,  or  insure 
that  there  are  no  highly  oscillatory  fast  modes,  but  these 
constraints  will  be  introduced  later  as  they  are  needed. 

Since  our  initial  purpose  is  only  to  identify  the  slowly 
varying  and  rapidly  varying  parts  of  system  (1.1),  the  two- 
time-scale  assumption  will  suffice. 

It  is  possible  for  the  spectrum  of  A to  cluster  by 
modulus  into  three  or  more  groups  so  the  system  dynamics  can 
be  analyzed  in  three  or  more  time  scales.  In  such  cases. 


' 
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the  system  possesses  two  or  more  small  parameters  defined 
analogously  to  (1.5).  Only  the  two-time-scale  case  will  be 
considered  here,  although  by  repeated  application  of  these 
time  scale  decoupling  techniques  the  multiple  time  scale  case 
could  be  treated. 

Given  a two-time-scale  control  system  with  n^  slow 
modes  and  n2  fast  modes,  we  partition  the  state  vector  x 
into  subvectors  x^  and  x2  of  dimensions  n^  and  n2# 
and  partition  the  A and  B matrices  accordingly  so  that 
the  control  system  (1.1)  can  be  rewritten  as 


t- 

t- 

A12_ 

"X1 

_Bll 

+ 

A21 

A22_ 

-2_ 

B2_ 

(1.6) 


Consider  a nonsingular  linear  transformation  of  the 
state  variables  x to  new  dynamic  variables  £ represented 
as 


£ = 


= rx. 


The  new  control  system  defined  by  this  transformation  of 


variables  is 


6 


l = 65  + 4>u  (1.7a) 

w = Y?  + Du  (1.7b) 

where  0 = TAr"1,  * = TB,  and  y = Cr”1.  If  system  (1.7) 
has  the  block-diagonal  form 

““  ®H^1  ^1^  (1.8a) 

£ 2 ~ 922^2  + ^2U  (1.8b) 

where  0 and  <p  are  partitioned  compatibly  with  (1.6), 
there  is  no  interaction  between  and  £2»  so  we  have 

achieved  a decoupling  of  the  system's  variables.  Since 
matrices  A and  0 are  similar,  they  have  the  same  eigen- 
values, i.e.. 


X (0)  = X (A)  . 

If  transformation  T decouples  system  (1.1)  and,  in 
addition,  the  eigenvalues  are  partitioned  >s 

X(0n)  = S,  X ( 022)  ° F'  (1.9) 


then  we  have  achieved  a time  scale  decoupling  of  the 
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\) 


I 

s 


) 


two-time-scale  system.  The  subsystems  (1.8a)  and  (1.8b)  will 
represent  the  system's  slow  dynamic  and  fast  dynamic  com- 
ponents, respectively.  Now  partition  r and  r”1  as 


r"1  = [ax  a21 


where  r^,  r2»  A-^  and  A2  have  dimensions  n^  x n, 

n2  x n,  n x n^,  and  n x Given  any  initial  state  x(0) 

and  control  function  u(t) , the  solution  to  the  initial 
value  problem  for  the  system  (1.1)  is  given  by 


x(t)  = eAtx(0)  + 


f. 


eA(t"T)B(T)u(T)dT. 


Employing  the  block-diagonal  structure  of  the  transformed 

system  (1.8),  the  solution  x(t)  can  be  defined  as  the 

direct  sum  of  slow  mode  and  fast  mode  constituents  x and 

s 

xf , i.e. , 


x(t)  = xg(t)  + xf(t) 


where 


0ii  t ft  0. . (t-T) 

xg(t)  5 A^e  J"L  rix  (0)  + j^e  11  I^B  (t)  u (t)  dx>  (1. 10) 


As  long  as  (1.8)  and  (1.9)  are  satisfied,  the  slow  mode  and 
fast  mode  constituents  xg(t)  and  xf(t)  are  unique  on 
[0,T]  for  all  r. 

Although  many  different  time  scale  decoupling  trans- 
formations T are  possible,  two  specific  choices  for 
which  have  especially  suitable  algebraic  structure  will  be 
considered  here.  These  two  transformations  are  presented  in 
Chapter  2. 


CHAPTER  2 


MODAL  DECOUPLING  AND  REDUCED  ORDER  MODELING 

One  well  established  method  of  time  scale  decoupling 
is  through  the  use  of  modal  variables  as  described,  for  ex- 
ample, by  Zadeh  and  Desoer  (1963)  or  Porter  and  Crossley 
(1972) . In  order  to  simplify  the  review  of  this  method,  as- 
sume for  the  present  that  A is  nondefective,  i.e.,  that  it 
can  be  put  into  diagonal  form  by  a similarity  transformation. 

2.1.  A Review  of  Modal  Variables 
Let  the  A matrix  be  represented  in  the  spectral 

form 

A = MJQ.  (1.1) 

Here  J is  a diagonal  matrix  whose  diagonal  elements  are  the 
eigenvalues  in  X (A) , M is  a modal  matrix  whose  columns  are 
the  corresponding  (right)  eigenvectors  of  A,  and  Q = M 
The  spectral  form  of  a general  matrix  A with  possibly 
generalized  eigenvectors  will  be  considered  in  Chapter  3. 

For  the  classical  method  of  modal  variables,  the 
transformation  T is  given  via  its  matrix  Q,  and  modal 
variables  are  defined  by 


9 


The  original  system  (less  output  variables) 


x = Ax  + Bu  (1.3) 

is  transformed  by  (1.2)  to  the  familiar  completely  decoupled 
form 


z = Jz  + QBu  (1.4) 

or,  equivalently,  the  n scalar  equations 

^ = Xj^z^  + q^Bu,  i = 1,  ...,  n 
T 

where  z = (z.^  ...,  z^)  , j = diag(X1,  ...,  Xn)  and  q.^ 
is  the  ith  row  of  matrix  Q.  This  completely  decoupled 
form  achieves  more  separation  of  eigenvalues  than  is  needed. 
That  is,  for  the  study  of  two- time -scale  systems  it  is  only 
necessary  to  separate  slow  and  fast  eigenvalues.  Modal 
variables  are  also  difficult  to  use  in  practice  for  the 
following  reasons: 

1)  The  physical  insight  which  may  be  inherent  in  the 
original  variables,  x,  may  be  lost  with  the  modal 
variables,  z. 


11 


2)  Modal  variables  may  be  complex  quantities  while  the 
physical  variables,  x,  are  most  apt  to  be  real. 

3)  The  transformation  Q is  nonunique,  so  the  modal 
variables,  z,  are  nonunique. 

4)  The  modal  matrix  M becomes  ill-conditioned  (cf.  Stewart, 
1973,  p.  297)  as  the  eigenvalues  of  A approach  each 
other. 

5)  The  accuracy  of  the  modal  variables  obtained  is  limited  by 
the  accuracy  of  the  numerical  scheme  used  to  calculate 
eigenvalues  and  eigenvectors. 

Nonetheless,  the  modal  form  (1.4)  is  sometimes  em- 
ployed in  practical  control  problems  since  it  may  yield 
useful  information  on  the  control  structure.  That  is,  the 
jth  element  of  the  row  vector  q^B  described  above  provides 
a measure  of  the  influence  of  the  jth  control  variable  upon 
the  ith  mode  of  the  system.  This  yields  insight  for 
controller  design. 

2.2.  The  LK  Decoupling  Transformation 
A new  transformation  useful  for  two-time-scale  systems 
was  introduced  by  Kokotovic  (1975)  and  provides  an  alterna- 
tive to  the  use  of  modal  variables.  This  transformation, 
which  will  be  referred  to  as  the  LK  transformation,  in- 
volves a two  step  reduction  of  the  system  matrix  A first 
to  block- triangular  form  and  then  to  block-diagonal  form.  An 


12 


earlier  form  of  this  transformation  appears  in  the  singular 
perturbations  literature,  cf.  Wasow  (1976),  O'Malley  (1969), 
and  Chang  (1972) . 

Consider  matrices  L and  K of  dimensions  n 2 x n.^ 
and  n^  x n 2 respectively  which  define  new  vectors  y^  and 
°f  dimensions  n^  and  , as  follows: 

Y1  = xx  + Ky  2 
y2  = Lx!  + x2 » 


or,  in  matrix  form 


13 


y 


Tx  where  T 


ri  + KL  fTl 


L I 


Here  I represents  the  identity  matrix  of  appropriate 
dimension  (n^  x n^  or  n2  x n^)  . 

If  the  L matrix  is  a solution  of  the  nonsymmetric 
algebraic  Riccati  equation  (ARE) 


I*H  ' A22L  - “l2L  + A21  - °<  (2-3) 

where  the  A^^  are  submatrices  of  the  A matrix  [as  de- 
fined in  Chapter  1,  equation  (1.6)  i.e.,  (1.1.6)]  then  the 
transformation  (2.1)  reduces  the  control  system  (1.3)  to 
block-triangular  form 


*1 

* X 

11  12 

xl 

_B1 

= 

+ 

j2_ 

1 o 

to 



_y2_ 

LBi  + B* 

where 


^11  E A11  ” A12L'  ^12  " A12  and  ^22  = A22  + ^12* 
Furthermore,  if  the  K matrix  satisfies  the  Lyapunov 


equation 


^22  — ^11^  ^12  — (2.4) 

then  the  transformation  (2.2)  decouples  (1.3)  yielding 

yx  = (2.5a) 

y2  = *22y2  + &2U'  (2‘5b> 

which  in  matrix  form  is  block-diagonal,  i.e.. 


y 


it? 

i 

i 

O 

I 

hi 

ii 

i 

y + 

o 

& 

22 

2 

Xy  + &u 


where  = (I  + KL)B1  + KB2  and  &2  = LB^  + B2. 

If,  in  addition  to  satisfying  ARE  the  L matrix 
separates  the  eigenvalues  as  required  by  (1.1.9),  that  is. 


Mftn)  = s,  x (k22)  = F 


then  the  block-diagonal  system  (2.5a,b)  is  time  scale 
decoupled.  Comparing  with  equation  (1.1.8),  we  find  that 
the  n..th  order  system  (2.5a)  represents  the  system's 


slow  mode  components  and  the  n2_t^  order  system  (2.5b) 
represents  its  fast  mode  part.  In  Chapter  3 we  will  show  that 
such  L and  K matrices  always  exist  for  weakly  two-time- 
scale systems  and  they  are  real  and  unique.  Accurate 
numerical  values  for  L and  K can  be  obtained  as  shown  in 
Chapter  4. 

One  attractive  feature  of  the  transformation  of  vari- 
ables defined  by  (2.2)  is  that  the  inverse  transformation  is 

x = 

Since  the  inverse  is  explicitly  determined  by  L and  K,  no 
numerical  matrix  inversion  procedure  need  be  used,  and  the 
inverse  exists  whenever  T is  defined.  Although  the 
physical  interpretation  of  the  modal  variables  z is  obscure, 
the  relationship  between  x and  the  decoupled  variables  y 
is  relatively  transparent.  If  we  decompose  the  state  vari- 
ables into  fast  and  slow  parts  xg  and  x^  as  defined  by 
(1.1.10),  we  will  obtain 


-K 


-L  I + LK 


y = 


T_1y- 


(2.6) 


*1 

xi 

Is 

xlf 

X = 

= 

+ 

_X2_ 

CO* 

CM 

*i 

*2t_ 

x 


(2.7) 
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and  since  x^  = y^  - Ky 2,  it's  apparent  that  x^g  = y^.  The 
variables  y^  therefore  have  the  same  physical  dimensions  as 
x^.  Nearly  the  same  analogy  can  be  made  between  y2  and 
x2,  i.e.,  x2  = -Ly1  + (I  + LK) y2  implies  x2f  ^ y2  when- 
ever LK  is  of  small  norm. 

To  summarize,  the  LK  transformation  has  several 
distinct  advantages  compared  to  the  use  of  modal  variables 
z and  the  modal  transformation  M described  in  the  previous 
section,  namely: 

(1)  The  physical  insight  inherent  in  the  x variables  is 


(2)  The  LK  transformation  T is  always  real  and  unique 
and  its  inverse  can  be  obtained  explicitly  in  terms  of 
L and  K without  using  numerical  matrix  inversion. 

(3)  Matrices  L and  K are  obtained  from  algebraic  equa- 
tions (2.3)  and  (2.4),  whereas  M and  Q have  to  be 
obtained  by  numerical  eigenanalysis  and  matrix  inversion. 

The  practicality  of  the  LK  transformation  will  be  fully 
demonstrated  in  later  chapters. 

Once  the  appropriate  L and  K matrices  have  been 
obtained  the  integrated  form  of  the  time  scale  decoupled 
system  (2.5)  becomes 

x ( t ) = T~1{e^tTx (0)  + fte^(t_x) 

'0 


Tb (t) u (t) dt } 


i 
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or,  applying  (2.7)  and  the  expanded  form  for  T: 


xs(t) 


(2.8) 


-L 


xf  (t) 


& t 

{e  [ (I  + KL)  K]  x (0)  + 


t (t-x) 


(x)  u (t)  dx } 


(2.9) 


-K 


I + LK 


K2t  ft  'K  (t-x) 

(e  [L  I]  x (0)  + j e &2  (x ) u (x ) dx  } , 


No  approximations  have  been  made  at  this  point  so  this  solu- 
tion x(t)  = x (t)  + xc(t)  is  exact, 
s f 

Without  additional  assumptions  on  x(0)  and 
B (t) u (t) , the  slow/fast  terminology  continues  to  be  merely 
suggestive  and  xg(t)  and  xf(t)  need  not  be  clearly 
distinguishable  as  slow  and  fast. 

The  decoupled  form  (2.8),  (2.9)  may  be  advantageous 
for  numerical  integration  since  the  slow  and  fast  modes  are 
integrated  separately  to  be  later  recombined  to  obtain  the 
state  x.  If  n^  ^ then  approximately  half  as  many 

multiplication  and  addition  operations  are  required  to 
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evaluate  derivatives  at  one  time  step  for  the  block-diagonal 
matrix  form  (2,5)  compared  to  the  original  (1.3). 

One  of  the  main  applications  of  modal  decoupling  is 
the  derivation  of  reduced  order  models  which  approximate  the 
slow  dynamic  behavior  of  the  original  system  (1.1).  Here 
maximum  benefit  is  achieved  if  n^  <<  n 2 so  that  the  reduced 
order  model  is  of  considerably  lower  order  than  the  original 
model. 


2,3.  The  Reduced  Order  Models 
of  Davison  and  Marshall 


In  order  to  be  able  to  neglect  fast  mode  behavior  and 

1 


have  x(t)  ^ x (t)  for  t >> 
= s 


we  will  need  the  fol- 


lowing fast  mode  stability  assumption: 


Is  I <<  -Re(f.)  for  all  f.  in  F.  (3.1) 

1 n^^ 1 j ] 

This  assumption  guarantees  that  the  fast  modes  are  all  stable 
and  well  damped  relative  to  the  slow  modes.  Since  it  is 
possible  for 


\sn^\  <<  | He (f j) | <<  |Im(f j) | , 

the  fast  modes  could  still  be  highly  oscillatory.  Under 
fast  mode  stability  if  we  are  to  consider  control  system 
response  for  large  times,  then  to  insure  bounded  state  x(t) 
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the  following  slow  mode  stability  assumption  will  also  be 
needed: 


Re(s^)  < 0 for  all  s^  in  S.  (3.2) 

The  two  conditions  (3.1)  and  (3.2)  taken  together  will  be 
referred  to  as  the  stability  assumption.  In  this  and  the 
following  section  it  is  assumed  that  A satisfies  the  fast 
mode  stability  assumption. 

Reduced  order  modeling  plays  a central  role  in  the 
study  of  large  scale  (high  order)  systems.  One  order  re- 
duction technique  based  on  modal  analysis  [suggested  by 
Nicholson  (1964)  and  presented  by  Davison  (1966)]  has  been 
widely  studied  (Sandell  et  al.  1978) . The  Davison  model  will 
be  derived  using  a notation  which  is  more  concise  than  the 
original  and  is  similar  to  that  of  Sandell  et  al.  (1978) . 
Consider  the  spectral  form  of  A (assumed  nondefective) 


A = MJQ 


partitioned  compatibly  with  system  (1.1.6)  as 


*11  A12 
_A21  A22_ 


11 

21 


(3.3) 
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where  A.^,  Mll'  J1  anc^  have  dimensions  n^  x n^, 

etc.  and 

A (Jj^)  = S,  X(J2)  = F.  (3.4) 

If  the  modal  variables  z are  also  partitioned  into  sub- 
vectors z^  and  z2  of  dimensions  n^  and  n2  then  the 
system  equations  in  modal  form  (1.4)  can  be  written  as  slow 
part 


Z1  = Jlzl  + 


[Q 


11 


Q12]Bu 


(3.5a) 


and  fast  part 


Z2  = J2Z2  + tQ21  Q221Bu*  (3.5b) 

The  transformation  from  modal  variables  back  to  state  vari- 
ables can  be  written  as 

X1  * Mllzl  + M12z2 

(3.6) 

x2  - M2l2l  + M22z2. 

The  Davison  model  is  based  on  the  approximation  that 
z2(t)  = 0 for  all  time,  i.e.,  the  fast  modes  can  be 


21 


disregarded.  This  implies  that  both  z2(t)  and  the  non- 
homogeneous  term  in  (3.5b)  are  small,  i.e.,  z2  ^ 0 and 
^21  ^22^ Bu  = This  may  or  may  not  be  a reasonable  ap- 

proximation depending  upon  the  particular  problem.  We  define 


X1  = Mllzl  ~ X1 


and 


(3.7) 


X2  * M2izi  = M21Mll_1^i  E x' 


-1 


assuming  exists.  Substituting  for  z in  (3.5a) 

yields  the  Davison  reduced  order  system 


-1^ 


X1  = M11J1M11  X1  + MlltQll  Q12]Bu(t) * (3.8) 


The  general  solution  of  (3.8)  can  be  written  as 


x(t)  = 


% 

x1(t) 


x2(t) 


(3.9) 


M 


11 


M 


21 


^ ^ rC  J . i 

M11  xi^°^  + J e ^Qn  Q12]B(T)u(T)dT} 


rt  J, (t-T) 
e 

0 
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which  is  Davison's  approximate  solution  for  the  original 
system  (1.3).  If  the  initial  state  x^(0)  is  chosen  as 

x^O)  = Q12]x(0) 

then  comparing  with  equation  (1.1.10),  it  is  apparent  that 

% 

x = x . However,  Davison  did  not  specify  this  value  for 

S 

% % 

x1(0),  but  assumed  x^tO)  = x^O). 

As  observed  by  Chidambara  (1967)  one  objection  to 
Davison's  model  (3.8)  is  that  for  the  case  of  constant  in- 
puts the  steady  state  values  of  x^  and  x2  are  both  in 
error.  That  is,  the  model  does  not  reproduce  the  steady 
state  values  of  the  original  system. 

Now  assume  that  A is  slow  mode  stable  as  well  as 
fast  mode  stable.  For  constant  input  u we  obtain  the 
steady  state  value  of  the  original  system  as 

x(»)  = -M 

In  contrast,  the  corresponding  steady  state  for  Davison's 
model  (3.8)  and  (3.9)  is  given  by 


Ji"1  0 


-1 


QB (») u 


(3.10) 
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^ i \ 

X(oo)  = - 

indicating  the  discrepancy.  (In  (3.10)  and  (3.11),  the  co- 
efficients are  evaluated  at  infinity.)  It  is  apparent  that 
(3.11)  is  simply  the  slow  mode  part  of  (3.10)  if  (3.10)  is 
written  in  expanded  form,  i.e., 

X (oo)  = Xg(oo)  . (3.12) 

A second  reduced  order  modeling  technique  developed 
independently  by  Marshall  (1966)  and  Chidambara  (1967)  is 
currently  being  used  to  reduce  the  order  of  large  scale 
systems  (cf.,  e.g.,  Skira  and  DeHoff  1977).  This  model  can 
be  derived  from  the  original  differential  equation  for  the 
partitioned  state  x^,  namely 

x.^  = + a12x2  + Biu*  (3.13) 

Using  the  inverse  modal  transformation  Q and  z = Qx,  we 
express  x2  in  terms  of  x^  and  z2  as 

X2  = “°22  1q21X1  + °22  1z2* 


[M 


11 


M 


21 


Jl_1(Qn  Q12]B(“)u 


(3.11) 
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The  invertability  of  Q22  follows  from  the  assumed  in- 
vertability  of  M^.  Now  the  approximation  z2  ^ 0 is  made 
in  (3.5b)  yielding 


z2(t) 


Q22]Bu(t) . 


Substituting  for  x2  and  z2  in  (3.13)  yields  Marshall's 
reduced  order  system 


X1  (A11  " A12Q22  Q21)x1 


+ [B1  - A12Q22  J2  [Q21  Q22JB]u(t), 


(3.14) 


with  x^(0)  = x^(0).  Marshall's  model  (3.14)  does  not  neces- 
sarily imply  that  J2Z2  anc*  ^21  Bu  are  smaH'  but 

only  that  there  is  no  difference  between  them.  This  is  the 
same  assumption  to  be  made  below  in  defining  the  zeroth  order 
outer  solution  of  a corresponding  singularly  perturbed  system. 

Writing  A^  and  A^2  in  their  spectral  forms 


MllJlQlj 


+ M 12J2Q2j ' 


j = 1,  2 


and  using  the  identity 


Q11  “ Q12Q22  1q 


21 


25 


it  follows  that 


(An  " A12Q22  Q21}  M11J1M11  ^ 


Comparing  (3.8)  and  (3.14)  shows  a strong  similarity  between 
Davison's  model  and  Marshall's,  i.e.,  they  have  the  same 
eigenvalues  and  eigenvectors.  The  coefficients  of  the  con- 
trol term  in  the  two  models  are  not  equal,  however,  i.e.. 


Mii[Qn 


°12>B 


* Bi  - a12°22'1j; 


-1 


[Q 


21 


Q22]B 


so  the  variable  x^  in  Marshall's  model  does  not  equal  x^g. 
The  primary  advantage  of  Marshall's  model  is  that  it  does 
yield  the  correct  steady  state  value 


(“)  = x1(«) 


(3.15) 


since 


x1(®) 


-MuJl*1Mil-1[Bi  - A12Q22*1J2‘1[Q21  Q22lBja, 


and  after  considerable  manipulation  it  can  be  shown  that 
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X1  (oo) 


[M11  M12] 


-1 


QBu 


which  agrees  with  (3.10).  (Coefficients  are  again  evaluated 
at  infinity.) 

Another  increasingly  used  approach  to  reduced  order 


1 

modeling  of  control  systems  is  the  reduced  order  approxima- 
tion of  singular  perturbation  theory.  In  this  method,  it 


is  assumed  that  a small  parameter  e has  already  been  identi- 
fied, and  the  control  system  can  be  written  in  the  form 


x^  = + A|2x2  + B|u  (3.16a) 

ex2  = A^Xj^  + A^2x2  + B2u.  (3.16b) 


Here,  the  prime  notation  is  used  to  indicate  that  a different 
partitioning  of  the  original  system  (1.3)  may  be  chosen.  The 
value  e must  be  small  relative  to  one  and  all  matrices  in 
(3.16)  are  usually  considered  to  be  the  same  order  of  magni- 
tude in  some  suitable  norm.  The  reduced  order  approximation 
to  (3.16)  obtained  by  setting  e = 0 and  assuming  A22 
invertible  is 


*1  “ ^Ail  “ A12A22  A21^X1  + (B1  ” A22  B2*U  (3.17a) 
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x2(t)  = -A^2~1(A^1x1(t)  + B£u(t)).  (3.17b)  l 

J 

In  order  for  the  model  (3.17)  to  yield  a good  approximation 
to  (3.16)  outside  boundary  layer  regions,  appropriate 
stability  hypotheses  are  required  (cf.  Kokotovic,  O'Malley 
and  Sannuti  1976) , like  the  fast-mode  stability  assumption  of 
the  previous  section.  One  of  the  key  advantages  of  this  ap- 

i 

proach  is  that,  if  the  system  can  be  written  in  the  form 

(3.16)  then,  even  for  systems  of  large  order,  the  reduced 
order  approximation  (3.17)  is  easily  obtained.  Equations  in 
form  (3.16)  are  usually  obtained  from  an  analytical  deriva- 
tion where  e may  represent  specific  physical  parameters  or 
ratios  of  physical  parameters,  though  such  scaling  is  usually 
for  low  order  systems.  In  some  studies  of  large  scale  problems 
where  the  system  equations  in  form  (1.3)  are  obtained  through 
numerical  methods,  the  partition  of  x into  x^  x2  is 
based  on  engineering  judgment  and  the  reduced  order  model 

(3.17)  is  obtained  without  considering  the  actual  value  of 
e (cf.  Kelley  1973;  and  others).  Harvey  and  Pope  (1976) 
employ  the  above  method  (termed  by  them  "residualization" ) 
to  reduce  a high  order  model  of  aircraft  longitudinal 
dynamics  for  the  purpose  of  evaluating  control  synthesis 
techniques . 


2.4.  A New  Reduced  Order  Model 
from  Singular  Perturbations 

The  singular  perturbation  method  can  be  readily 

applied  to  the  two-time-scale  control  system  in  block- 

diagonal  form,  namely 


^1  = ^llyl  + ^lu  (4.1a) 

y2  ' ^22y2  + &2U*  (4'lb) 


In  order  to  clearly  describe  the  singular  perturbation  nature 
of  (4.1),  introduce  the  new  small  parameter 


Given  any  small  number  6 > 0,  the  norm  |]  • ||  can  be  chosen 
so  that 


«*11«  i KJ  + Slfn2' 

and  for  any  compatable  norm 


(cf.  Stewart  1973,  p.  184,  270)  so  that  a 
be  smaller  than  y,  cf  (1.1.5),  i.e.. 


will  generally 
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o < p + 6 < y + 6 <<  1. 

n2 

Multiplying  (4.1b)  through  by  either  small  parameter  a or 
y clearly  shows  its  singularly  perturbed  nature.  Indeed, 
if  one  is  to  consider  a linear  control  system  as  a singular 
perturbations  problem,  then  with  no  loss  of  generality  the 
system  may  be  considered  in  the  form 

Y i = + &xu  (4.2a) 


^2  = *22^2  + V 


(4.2b) 


where  A22 
value  of 


= P^22  and 

'X,  'Vi 

y,  A or  B. 


B2  e y&2  with  no  ambiguity  as  to  the 
Note  that  this  rescaling  makes  the 


entries  in  A22  of  order  0(1). 

The  singular  perturbation  procedure  eliminates  the 
derivative  term  in  the  limit  U -*■  0 because  that  term  is 
multiplied  by  the  specific  small  parameter  y.  In  its 
physical  interpretation,  the  derivative  term  is  omitted  be- 
cause the  dynamics  of  are  fast  and  the  derivative 

y2(t)  is  small  after  a typical  fast  dynamic  transient  has 
died  out.  Equation  (4.2)  illustrates  the  equivalence  of 
these  two  approaches. 

The  following  new  reduced  order  model  approximating 
the  system's  slow  dynamics  is  obtained  by -setting  y = 0 


| 

) 
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in  (4 . 2b) : 


y-L  = (4.3a) 

y2  E ~^22_1^2U  (4.3b) 


x(t) 


-K 

yx(tr 

I + LK 

y2  (t) 

(4.3c) 


(The  last  notation  corresponds  to  the  new  reduced  order 
model.)  For  initial  value  problems,  the  initial  conditions 
appropriate  for  (4.3a)  are 

y^O)  = (I  + KLlx^O)  + Kx2(0).  (4.3d) 

(If  our  system  were  fast  mode  unstable,  we  could  analogously 
specify  terminal  values  for  the  state  vector  and  require  that 
the  control  component  B(t)u(t)  on  [0,T]  not  excite  these 
unstable  modes.  The  more  difficult  problem  when  ^22 
conditionally  stable  deserves  further  study  with  both 
boundary  values  x(0)  and  x(T)  and  control  term  B(t)u(t) 
appropriately  restricted  [O'Malley  1969;  O'Malley  and 
Anderson  1979].)  Comparing  this  model  with  the  original 
decomposition  of  x into  slow  and  fast  components 


r 

1 
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x = xg  + Xf  we  have  an  exact  representation  of  the  slow 
component , i . e . , 


1 

xs(t)  = 

-L 

t-* 

rt 

II 

-L 

{e 


&ut 


[ (I  + KDx^O)  + Kx2  (0)  ] 


+ 


,t  n (t-x) 
e 11  [I  + KL 

0 


K]B(x)u(x)dT}. 


(4.4) 


Assuming  the  system  (1.3)  is  also  slow  mode  stable, 
the  steady  state  response  of  (4.3)  for  constant  input  u 
is 


x(°°)  = 


I -K 


-L  I + LK 


"-ft  -i 
11 


-A  -1 
22 


I + KL  K 


Bu 


= -T"1A~1TBU 


(4.5) 


(with  evaluation  at  infinity) . Comparing  this  form  with  the 
steady  state  response  of  (1.3),  i.e.. 


x(<*>) 


A 

B (")  u 


/ 


the  equivalence  of  the  steady  state  responses 
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X (°°)  = X(oo) 

is  shown  by  the  relationship  between  the  original  system  and 
the  block  diagonal  form,  i.e.. 


A = T 


Thus  the  new  reduced  order  model  (4.3)  reproduces  the  slow 
mode  component  of  the  system  response  to  both  initial  condi- 
tions and  control  inputs,  and  yields  the  correct  steady  state 
response  to  constant  inputs  overcoming  both  difficulties  en- 
countered in  the  Davison  and  Marshall  models. 

If  the  original  control  system  includes  output 
variables  w,  i.e.. 


x = Ax  + Bu 


w = Cx  + Du, 


(4.6) 


then  the  new  reduced  order  model  may  be  written  in  a similar 
compact  form.  If  a new  output  matrix  & is  defined  as 


2*  = C T-1 


and  partitioned  as 
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then  the  new  model  (4.3)  may  be  written  as 


Yl  = + V 

W = + fcu 


(4.7) 


where  fr  = D - (^^22  ^2  initial  condition  y^(0) 

specified  by  (4.3d). 


2.5.  A Second-Order  Example 
The  differences  in  the  accuracy  of  approximation  be- 
tween Davison's  model,  Marshall's  model  and  the  new  model 
(4.3)  can  be  illustrated  by  the  following  simple  second- 
order  example  used  by  Chidambara  (1967)  and  Wilson,  Fisher 
and  Seborg  (1972)  to  test  order  reduction  methods: 


0 1 


-4  -5 


u. 


(5.1) 


For  Davison's  method  (3.8)  the  reduced  order  model  is 


yx  = -yL  + u/3 

y1(0)  = 4x1(0)/3  + x2(0)/3 

y2 (t)  = u ( t) /4 ■ (5.4) 

xl(t)  = yx(t)  - y2(t)/3 
x2(t)  = -y1(t)  + 4y2(t)/3 

The  response  of  the  three  models  to  initial  condi- 
tions (x1(0) ,x2 (0) ) = (1,1)  is  illustrated  in  Figure  2.1. 
The  new  model  is  clearly  superior  here  as  (4.3d)  provides 
the  true  projection  of  initial  conditions  x(0)  onto  the 


I = True  solution 


II 

III 

IV 


Davison' s 
Marshall ' 
New  model 


model 
t model 


I & III 


slow  eigenspace  (i.e.,  the  dimensional  space  spanned  by 

the  slow  eigenvectors)  whereas  neither  of  the  other  models 
provide  this  feature.  Since  the  LK  transformation  T is 
a time-scale-decoupling  transformation  r (cf.  (1.1.7)) 
which  exactly  decomposes  n-dimensional  space  into  fast  and 
slow  eigenspaces,  the  new  model  correctly  projects  the  slow 
mode  y1(t)  back  onto  the  original  state  x(t). 

The  responses  of  the  three  models  to  unit  step  input 
u(t)  = 1,  ramp  input  u(t)  = t,  and  sinusoidal  input 

TT 

u(t)  = sin(jt)  are  illustrated  in  Figures  2.2  through  2.4. 

The  boundary  layer  structure  of  the  fast  mode  corresponding 
to  eigenvalue  -4  can  be  clearly  seen  in  Figures  2.1  and 
2.2  as  the  difference  between  the  true  solution  and  the  new 
model.  In  the  first  three  cases  the  new  model  provides  a 
more  accurate  approximation  to  the  true  solution  than  either 
Davison's  or  Marshall's  models.  For  the  case  of  sinusoidal 
input  the  new  model  and  Marshall's  model  lead  and  lag  the 
true  response,  respectively,  and  might  be  considered  equally 
accurate.  Note  that  in  all  cases  since  the  new  model  approxi- 
mates the  fast  mode  as  asymptotically  fast,  i.e.,  the  time 
scale  in  (4.2b)  shrinks  to  zero,  the  new  model  leads  the  true 
response.  Since  these  models  are  only  intended  to  represent 
the  slow  mode  dynamics  of  the  system,  they  can  only  be 
expected  to  be  accurate  for  low  frequency  inputs,  i.e., 
frequencies  less  than  the  order  of  magnitude  of  the  time 


varialbe 
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I = True  solution 
II  = Davison's  model 
III  = Marshall's  model 


Figure  2.4.  Response  of  reduced-order  models  to  a sinusoidal 
input  for  Example  2.5. 
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constant  1/ | s |.  The  new  order  reduction  method  and 
nl 

Marshall's  method  will  be  compared  again  in  Chapter  5 when 
they  are  employed  to  model  the  slow  mode  dynamics  of  an  F8 
aircraft. 


2.6.  Feedback  Matrix  Design 
via  the  LK  Transformation 

Often  dynamic  systems  are  modeled  as  linear  control 
systems  for  the  purpose  of  designing  a feedback  controller  to 
achieve  specific  closed-loop  eigenvalue  locations.  As  will 
be  shown,  this  task  can  be  greatly  simplified  by  use  of  a 
time  scale  decoupling  transformation  such  as  the  LK  trans- 
formation. 

It  is  our  purpose  here  to  show  that  given  any  time- 
scale  decoupling  transformation  r,  such  as  T,  an  n- 
dimensional  eigenvalue  placement  problem  can  be  reduced  to 
separate  eigenvalue  placement  problems  of  dimensions  n^ 
and  n2»  Here  the  control  matrix  B and  the  output  matrix 
C will  be  taken  to  be  constant,  and  n^  and  n2  could  cor- 
respond to  any  weakly  two-time-scale  division  of  eigenvalues. 

It  is  well  known  that  if  a constant  linear  control 
system  is  controllable,  i.e., 

rank  (B  SAB|  . ..  JAn  ^B)  = n,  (6.1) 

then  there  exists  at  least  one  real  m * n dimensional 
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feedback  matrix  H such  that  the  closed-loop  eigenvalues 
given  by 

X(A  + BH)  (6.2) 

can  be  placed  arbitrarily  as  long  as  complex  eigenvalues 
appear  in  conjugate  pairs. 

If  the  original  system  is  controllable  then  the 
decoupled  slow  and  fast  systems  (2.5)  are  also  controllable. 
This  can  be  demonstrated  as  follows.  Let  the  block- 
decoupling transformation  T premultiply  the  matrix  (6.1) 
to  yield 


[tb;  (tat'1) tb;  ...  ;(tat-1)  ...  (tat'1)tb) 


If  the  n x nm  matrix  (6.1)  is  of  full  rank,  then  the 


matrices 
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si! 

n_1 

...  ,An 

2^1  , 

(6.3a) 

'62i*22 

&2i 

*•*  J 22 

62] 

(6.3b) 

are  each  of  full  rank. 

and 

controllability  of 

the  pair  (A,B) 

implies  controllability  of  and  ^22'^2^’  However, 

the  converse  is  not  true  since  for  arbitrary  ^22' 

full  rank  matrices  (6.3a)  and  (6.3b)  could  have  linearly 
dependent  rows.  Conditions  under  which  controllability  of 
the  fast  and  slow  subsystems  implies  controllability  of  the 
original  system  are  described  by  Chow  (1977)  and  Sannuti 
(1978).  (If  the  pair  (At,Ct)  is  controllable  where  C is 
the  system  output  matrix  in  (4.6)  then  the  control  system 
(4.6)  is  observable.  If  the  original  control  system  is 
observable,  then  it  can  be  shown  by  duality  that  the  de- 
coupled control  systems  (4.1)  are  also  observable.)  For  the 
remainder  of  this  section,  we  will  assume  that  the  original 
system  (4.6)  is  controllable. 

The  design  of  a feedback  matrix  H to  achieve  a spe- 
cific closed-loop  eigenvalue  location  is  very  straightforward. 
Suppose  that  it  is  desired  to  relocate  the  n^  slow  eigen- 
values S to  n^  new  eigenvalue  locations  S'.  If  one  can 
find  an  m x n.^  feedback  matrix  which  satisfies 

M^ll  + = S' 


(6.4) 


then  the  system  (4.1)  with  feedback 
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implies  that 


(6.5) 


Thus  the  fast  eigenvalues  F are  unchanged  and  the  cor- 
responding feedback  matrix  H for  the  original  system  (4.6) 
can  be  found  as 


u = Hx  for  H = [2^  0]  T (6.6) 

since  y = Tx. 

If  in  addition  to  relocating  the  slow  mode  eigen- 
values S it  is  also  necessary  to  relocate  the  open-loop 
fast  eigenvalues  F to  closed- loop  values  F' , then  this 
can  be  done  as  follows.  After  solving  for  feedback  matrix 
in  (6.4),  find  the  unique  (assuming  S'  and  F have  no 
common  eigenvalues,  Gantmacher  1959,  p.  225)  n2  x n^  matrix 
P satisfying  the  Lyapunov  equation 
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r 

I 


f 

I 


i 


p(^H  + S^)  - &22p  + ^2^1  = 0>  (6.7) 

The  block-triangular  form  (6.5)  can  then  be  transformed  to  a 
new  block-diagonal  form  by  introducing  new  fast  variables 

Y 2 = Pyi  + y2'  (6,8) 

Now  relocate  the  fast  eigenvalues  F to  new  values  F'  by 
finding  an  m x feedback  matrix  fr ^ satisfying 

\(%22  + (P^x  + %2)&2)  = F'*  (6.9) 

The  appropriate  gain  matrix  for  the  original  state  variable 
x is  then  given  by 


u 


Hx 


(6.10) 


where  H = [ (ft^  - frjP)  T^us  the  block-decoupling 

transformation  T can  be  used  to  exactly  relocate  both  slow 
and  fast  eigenvalues  separately  via  state  feedback.  Ap- 
proximate methods  of  locating  slow  eigenvalues  (Georgakis  and 
Bauer,  1978) , or  separately  locating  slow  and  fast  eigen- 
values (Chow  and  Kokotovic,  1976)  based  on  the  singularly 
perturbed  system  (3.16)  have  also  been  proposed. 


t 

J 


x 

p 
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Many  variations  of  this  approach  are  possible  to 
facilitate  the  feedback  matrix  design  task.  For  example,  by 
repeated  application  of  the  block-decoupling  transformation, 
either  of  the  decoupled  subsystems  in  (4.1)  could  be  further 
transformed  into  block-diagonal  form  and  a complete  eigen- 
value relocation  could  be  achieved  in  three  consecutive  de- 
sign steps  analogous  to  the  above  two  design  steps  (6.4)  and 
(6.9).  In  many  practical  applications  (e.g.,  Skira  and 
DeHoff,  1977)  the  fast  modes  are  all  stable  and  well  damped 
as  defined  in  (3.1)  and  (3.2)  so  only  the  slow  mode  eigen- 
values need  be  relocated. 

The  practicality  of  this  procedure  depends  directly 
upon  being  able  to  calculate  the  appropriate  L and  K 
matrices.  As  will  be  shown  in  Chapter  4,  accurate  approxima- 
tions to  the  matrices  may  be  readily  computed. 


1 


CHAPTER  3 


EXISTENCE  AND  UNIQUENESS  OF  THE  L AND  K 
TIME-SCALE  DECOUPLING  MATRICES 

In  this  chapter  the  general  solution  to  the  non- 
symmetric  algebraic  Riccati  equation  (ARE),  (II. 2. 3)  is 
described.  As  shown  in  Chapter  2,  the  task  of  finding  an 
(n  - n^)  * n^  L matrix  such  that  transformation 


reduces  a partitioned  linear  system 


to  block-triangular  form 


A11  “ A12L  A12 


A22  + “12 
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is  equivalent  to  solving  the  ARE.  Past  work  on  the  existence 
of  solutions  to  ARE  is  reviewed,  and  general  conditions  for 
existence  of  L which  must  be  satisfied  by  matrix  A and 
by  the  integer  n^  (where  1 <_  n^  < n)  are  presented  in 
Theorem  3.1.  It  is  shown  that  L can  be  represented  in 
terms  of  eigenvectors  and  generalized  eigenvectors  of  A, 
and  that  there  is  a simple  relationship  between  the  eigen- 

'V/ 

values  of  A^  = A^  - A^L  and  these  eigenvectors. 

In  order  to  describe  the  solution  L to  ARE  we  first 
need  to  review  the  eigenstructure  of  square  matrices.  After 
the  solution  to  ARE  is  presented,  it  is  shown  that  the  matrix 
K satisfying  the  Lyapunov  equation  (II. 2. 4)  exists  whenever 
L exists.  In  Section  3.3  conditions  for  L and  K to  be 
real  are  described.  It  is  then  shown  in  Theorem  3.2  that 
for  a given  (weakly)  two-time-scale  system,  the  L matrix 
is  unique.  From  these  results  it  follows  that  time-scale- 
decoupling transformation  T (II. 2. 2)  exists,  is  real  and 
unique  for  any  given  (weakly)  two- time-scale  system. 

3.1.  A Review  of  Matrix  Eigenstructure 
At  this  point  it  is  not  assumed  that  A has  any 
special  properties  such  as  the  weakly  (or  strongly)  two- 
time-scale  property,  but  only  that  it  is  a real  square  matrix 
of  dimension  n x n with  p distinct  eigenvalues 
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For  all  matrices  A there  exists  a nonsingular  modal 
matrix,  M,  which  transforms  A to  Jordan  (canonical)  form 


M 1AM 


(1.1) 


(note  J.  is  defined  in  (II. 3. 3)  and  the  J.  are  m.  x m. 

1 ill 

Jordan  blocks  (Noble  and  Daniel,  1977,  p.  368)  . In  general, 
p k,  and  the  Jordan  block  J.  corresponding  to  the  eigen- 
value X.  will  be  indicated  as  J.,  . a Jordan  block  is  a 

3 i A. 

square  matrix  of  the  form 


1 

X . 

U 


(1.2) 


(If  nu  - 1,  the  one's  do  not  appear.)  The  number  k of 
Jordan  blocks,  the  numbers  n^,  ...,  m^  and  the  values 
Ar  ...»  Ap  are  uniquely  determined  by  A.  However, 
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the  ordering  of  Jordan  blocks  within  J is  arbitrary. 

If  A j is  an  eigenvalue  with  multiplicity  £ j , let 

J represent  the  l . x l . matrix  formed  from  the  Jordan 
J 3 3 

blocks  in  J with  eigenvalue  A . , and  let  M represent 

3 Aj 

the  n x matrix  formed  from  the  £j  columns  of  M cor- 

responding to  J.  , so  that 


AM.  = M.  J.  . (1.3) 

Aj  Aj  Aj 

The  -dimensional  eigenspace  of  Cn  spanned  by  the  columns 
of  is  uniquely  determined  by  A. 

j 

Taking  a complex  conjugate  (indicated  by  an  overbar) 
of  (1.3), 


= ma  ja 

Aj  Aj  Aj 


(1.4) 


implies  that  if  J.  , is  in  j,  then  J.  , is  also  in 

X'xj 

J.  Also,  without  loss  of  generality,  My  may  be  repre- 


sented as  M^  , assuming  the  ordering  of  blocks  within  J^ 

j j 

and  is  the  same. 

j 

Corresponding  to  each  Jordan  block  J.  . there  is  a 

l,Aj 

(right)  eigenvector  v,  . which  is  a column  of  M 

i'  3 


I 
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satisfying 


(A  - A . I)  v . . = 0. 

D Aj 


(1.5) 


If  A.  is  a distinct  eigenvalue,  then  v.  is  unique  up 

I 1 / A j 

to  a multiplicative  complex  scalar  constant.  This  non- 
uniqueness can  be  removed  by  normalization. 

If  Aj  has  multiplicity  greater  than  one,  then  any 
linear  combination  of  the  eigenvectors  of  Aj  is  also  an 
eigenvector.  If  iru  > 1,  then  corresponding  to  the  Jordan 
block  J.  there  are  generalized  (right)  eigenvectors 

X'Aj 

v.  . of  rank  r for  r = 2,  ...,  m.  which  satisfy 
1 , A j 1 


(A  - y >rvI,»j 


(a  - A,i)r"V:  f o, 

J 1 • A j 


(1.6) 


These  generalized  eigenvectors  are  also  columns  of  M cor- 
responding to  J.  . . Any  linear  combination  of  generalized 

i,Aj 

eigenvectors  of  lower  rank  may  be  added  to  a generalized 
eigenvector  and  (1.6)  will  still  be  satisfied.  One  straight- 
forward procedure  for  computing  the  generalized  eigenvectors 
is  described  by  Chen  (1970) . 
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3.2.  The  General  Solution  to  the 
Nonsyirunetric  Algebraic 
Riccati  Equation 

Given  matrix  A and  an  integer  n^  satisfying 
1 <_  n^  < n,  we  can  determine  the  transformed  matrix  ft  by 


ft  = 


ft  ft 
11  12 


ft  ft 
21  22 


I 0 


L I 


A11  A12 


A21  A22 


I 0 


-L  I 


(2.1) 


and  ft^  = 0 is  ecJu;*-va^ent  to  ARE 


^11  - A22L  - LA12L  + A21  = °' 


(2.2) 


In  his  previous  study  of  ARE,  Kokotovic  (1975)  used  a 'fixed 
point  theorem  to  show  that  if  the  norm  inequality 


l!A22_1|l  £T(HA11  - A12A22"1a2iII  + IIA12IIIIA22~1a2iH)~1 


is  satisfied  (assuming  A^2  exists)  then  there  exists  a 
real  solution  to  ARE  which  is  the  unique  limit  of  the 
sequence  of  matrices 


L0  A22  A2l' 


-1 


(2.3) 


Li+1  = A22  (A21  + LiAll  _ LiA12Li) ' 1 “ 0f  X'  2>  " 


i 

I 


i 
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In  another  study  Narasimhamurthi  and  Wu  (1977)  proved  that 
if  all  the  eigenvalues  of  A are  distinct,  then  there 
exists  a complex  matrix  satisfying  ARE.  Also,  Veljko  (1977) 
presented  the  solution  to  the  ARE  for  the  special  case  where 
each  multiple  eigenvalue  has  only  a single  eigenvector.  Both 
of  these  works  described  the  relationship  between  the  matrix 
L satisfying  ARE  and  the  eigenvectors  of  A.  These  results 
are  included  in  Theorem  3.1  which  describes  general  necessary 
and  sufficient  conditions  for  the  existence  of  a matrix  L 
satisfying  ARE.  The  key  results  of  Sections  3. 2-3. 4 were 
presented  in  Anderson  (1978) . 

Potter  (1966) , Martensson  (1971)  and  Kucera  (1972) 

have  also  described  the  eigenstructure  of  L for  the  special 

case  of  the  symmetric  algebraic  Riccati  equation  from  the 

optimal  linear  regulator  problem.  For  the  symmetric  Riccati 
. T 

equation  n^  = n^,  A^  = A^  and  Aj2'  A21  are  syiranetric 
and  positive  semi-definite.  In  the  large,  our  analysis 
parallels  theirs. 


Theorem  3.1.  Matrix  L satisfies  the  ARE  if  and  only  if 
(1)  some  subset  of  the  numbers  {nu  } can  be  chosen  so  that 


nu  ^ ...  m j n ^ • 


(2.4) 


Then 
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(2.5) 


is  an  x matrix  and 


Pi 


(2.6) 


is  an  x ^ matrix  where  i',  ...,  j'  are  the 
indices  between  1 and  k not  used  in  (2.5).  Also  the 
matrices  M and  Q = M ^ are  partitioned  compatibly 
with  and  3^  so  that 


A = mS'Q  = 


11 


21 


(2) 

(3) 


M.^  is  of  full  rank. 


“ll  = -M21* 


(2.7) 


(2.8) 


Proof.  Assume  first  that  L satisfies  ARE  (2.2).  There 
exists  an  n^  x n^  matrix  in  Jordan  form  and  a non- 

singular matrix  X such  that 
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A11  " A12L  = X6!3*”1  (2.9) 

and  ARE  can  be  rewritten  as 

A21  " A22L  = _L(A11  ~ A12L)  = “LXGjX”1*  (2.10) 

Postmult.iply  (2.9)  and  (2.10)  by  X and  let  Y = -LX  so 
that 


Thus  is  an  n^  x n^  block  of  the  Jordan  canonical  form 

of  A,  X is  of  full  rank,  and  the  L matrix  satisfies  the 
linear  equation  LX  = -Y  formed  from  the  appropriately  par- 
titioned columns  of  a modal  matrix  corresponding  to  G^.  As 
described  earlier,  the  modal  matrix  is  determined  only  up  to 
the  choice  of  basis  vectors  of  the  eigenspaces,  so  without 
loss  of  generality  we  can  choose  M^,  M2i  as  X'  This 

completes  the  first  half  of  the  proof. 

Now  assume  (1) , (2)  and  (3)  above  are  satisfied.  We 
need  to  show  that  the  lower  left  block  of  (2.1)  is 

zero.  Since  QM  = I,  it  follows  that 
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Q21M11  + Q22M21 


0 


so 


°21  “Q22M21M11 


-1 


and 


Q22L  °21‘ 


(2.11) 


Since  the  n2  rows  of  [Q21  Q22J  are  of  full  rank,  it  fol- 
lows that  Q22[L  ^ and  ^22  are  ^oth  °f  full  rank.  Write 
A as  (2.7)  and  complete  the  product  (2.1)  so  that 


M M 
11  n12 


M21  M22 


»1  ° 


o 3r„ 


Q11  °12 


°21  °22 


"°22  °21  1 


“llVll'1 

MllVl2  + M123'2Q22 

ft 

A11 

A12_ 

0 

°22  3f2Q22 

0 

W3 

to 



(2.12) 


since 


M11  “ Q11  “ Q12Q22  1q21 


(2.13a) 
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and 


M M 
21  11 


This  completes  the  proof. 


(2.13b) 


This  theorem  provides  the  key  result  which  will  be 
used  later  when  computing  the  L matrix;  namely,  if  L 
satisfies  (2.8)  then 


^11M11  MlA' 


(2.14a) 


and 


X Q -1 
22U22 


°22  ^2 


(2.14b) 


so  the  eigenvalues  and  eigenvectors  of  3^  and  are 

determined  via  3^,  and  3^,  Q^2  \ respectively. 

The  case  of  generalized  eigenvectors  causes  no  difficulty 
since  any  Jordan  block  appearing  in  also  appears  in  3^ 

and  conversely. 

Theorem  3.1  shows  that  under  the  very  special  condi-, 
tion  that  either  (1)  or  (2)  above  is  not  satisfied,  solutions 
to  ARE  may  not  exist.  However,  the  (weakly)  two- time-scale 
property  guarantees  that  (1)  can  be  satisfied,  and  if  one  is 


allowed  to  reindex  the  state  variables  x,  then  (2)  can 


also  be  satisfied. 

Recall  that  if  L satisfies  ARE  and  there  exists  a 
matrix  K satisfying 

^22  " ^11K  + A12  = 0 (2.15) 

then  TaT  1 = ^ is  block-diagonal.  Such  a K satisfying 
(2.15)  can  be  explicitly  represented  in  terms  of  the  modal 
matrix.  Since  MQ  = I, 


M11Q12  + 


M12Q22 


0, 


expressing  ^3/  ^22  anc^  A12  '*'n  sPectral  form,  it  follows 

that  (2.15)  is  satisfied  by 


K 


M11Q12 


12w22 


(2.16) 


i • 6*  1 


(_M12Q22)  (Q22  ^2°22)  " {M11^1M11  ^ (M11Q12) 


+ (m11^1Q12  + M12'^2Q22)  = °* 


Thus  K exists  whenever  L exists 


We  then  have  the  alternative  expressions 


L ‘ Q22"1«21  * -M21M11~1- 


(2.17) 


Applying  identities  (2.14)  we  can  express  the  LK  trans- 
formation T and  its  inverse  in  terms  of  partitions  of  the 
modal  matrix,  viz., 


+ KL  K MllQll  MUQ12 


1 °22  °21  1 


(2.18) 


r1  = 


I -K 


M12Q22 


-L  I + LK  M21Mu  M22Q 


. (2.19) 


3.3.  Conditions  for  Matrices 
L and  K to  be  Real 


One  of  the  primary  limitations  to  the  usefulness  of 
the  modal  matrix  M and  its  inverse  Q is  that  they  are 
complex  whenever  the  system  has  complex  eigenvalues.  As 
shown  below,  under  certain  restrictions,  the  L matrix 
satisfying  ARE  will  be  real. 
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In  order  to  facilitate  the  presentation,  we  will  as- 
sume that  the  Jordan  blocks  within  the  Jordan  form  3^  of 
(2.5)  are  ordered  so  that 


% 


J 

o 


J 


where  Jq,  J+  and  J_  correspond  respectively  to  Jordan 

blocks  in  3^  with  Im(Xi)  = 0,  Im(Ai)  > 0,  and  ImjA^ 

< 0,  and  the  dimensions  of  J , J,  and  J are  m x m , 

o + - o o' 

m+  * m+  and  m_  * m_,  with  mQ  + n+  + m = n^.  Also,  let 
the  mQ  columns  of  M corresponding  to  Jq  be  represented 
as  M and  likewise  for  M,  and  M so  that 

J 

o 

A [M  M,  M ] = [M  M M ] 
o + — o + - 


Then  if 


(3.1) 


the  L and  K matrices  satisfying  ARE  and  (2.15)  are  real 
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This  can  be  shown  as  follows.  Let  M_  be  repre- 
sented as  M+.  By  Theorem  3.1,  L satisfies  LMi;l  = -M21. 
Postmultiply  [HQ  M+  M+]  by  the  nonsingular  matrix 


where  for  clarity  the  dimensions  of  the  identity  matrices  are 
explicitly  noted,  to  obtain  the  real  matrix 


raL, 

M 

11 

_ 

11 

to 

1^ 

M 

21 

Thus  L satisfies 


1 


“11=1 


-M21=l 


or  = ”^2]  anc*  L therefore  real. 

If  the  Jordan  blocks  of  3^  appear  in  complex  con- 
jugate pairs,  cf.  (3.1),  then  so  do  those  of  By  the 

above  argument,  there  also  exists  a nonsingular  complex 
. ^ 

matrix  C2  of  dimension  n2  x n2  such  that 


K. 
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'Vi  ““ 

M12 

M 

12 

'Vi 

M22j 

M 

22 

is  real.  Equations  (2.16)  and  (3.2)  imply  that 


(3.2) 


K = -M12Q22  - M12(M22  - 


m12&2&2  (H22  - m12* 


-1..  ."I 


-Jl  (f)(  - ft  & ^ft  ) 

12lM22  21  11  W12 


so  K is  also  real  under  conditions  (3.1). 

From  these  results  it  follows  that  if  the  LK  trans- 
formation T decouples  the  time  scales  of  a (weakly)  two- 
time  scale  system,  then  T is  real. 

3.4.  Uniqueness  of  the  L 
and  K Matrices 

As  shown  by  Narasimhamurthi  and  Wu  (1977)  the  ARE 
will  have  a unique  solution  only  under  the  very  restrictive 
condition  that  the  quadratic  term  is  zero  and  the  equation 
becomes  linear.  In  fact,  we  will  show  below  that  there  are 
in  general  many  solutions  of  ARE.  However,  the  L and  K 
matrices  in  T which  decouple  a given  (weakly)  two-time-scale 


-1 
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system  so  that  slow  and  fast  modes  are  separated  are  unique. 
This  is  proven  in  the  following  theorem. 

Theorem  3.2.  For  a given  (weakly)  two-time-scale  system 
(with  n^  specified) , there  exists  at  most  one  decoupling 
matrix  L satisfying  both  ARE  and 

A (#2^)  = S,  X(&22)  = F.  (4.1) 

Proof . Assume  that  L and  L'  both  satisfy  ARE  and  (4.1). 
Applying  Theorem  3.1  there  exist  full  rank  matrices 


X 

X' 

and 

Y 

Y' 

such  that  LX  = -Y  and  L'X1  = -Y'.  Moreover,  the  columns 
of  (4.2)  are  (possibly  generalized)  eigenvectors  of  A cor- 
responding to  eigenvalues  S.  Since  there  is  a unique  n^ 
dimensional  eigenspace  spanned  by  these  (generalized)  eigen- 
vectors, there  exists  a nonsingular  matrix  u such  that 
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Postmultiply  L'X'  = -Y'  by  £ to  obtain  L'X  = -Y.  From 
Theorem  1,  X is  of  full  rank,  so  L = L'  = -YX-1,  com- 
pleting the  proof. 

The  uniqueness  of  the  K matrix  is  insured  by  the 
two-time-scale  property.  As  shown  by  Gantmacher  (1959) , the 
Lyapunov  equation  (2.15)  has  a unique  solution  provided 
and  ^22  have  no  common  eigenvalues. 

Since  L and  K are  uniquely  determined,  the  de- 
coupling transformation  T is  thereby  unique. 

Although  the  transformation  T is  unique,  there  are 
generally  many  real  solutions  to  ARE.  Thus  any  algorithm 
to  obtain  L should  avoid  converging  to  these  other  solu- 
tions. In  a special  case,  the  number  of  distinct  solutions 
can  be  enumerated. 

Corollary  3.3.  If  A has  (real)  distinct  eigenvalues  and 
each  set  of  n^  columns  taken  from  the  n columns 
[Mu  M12^  '*'s  rank>  then  there  are  (R^)  distinct 

(real)  solutions  of  ARE  corresponding  to  the  number  of  dis- 
tinct ways  of  partitioning  the  n eigenvalues  into  sets  of 
n^  and  n2  members. 

Proof.  By  Theorem  3.1  and  3.2  there  exists  one  solution 
for  each  of  the  (n^)  partitions  of  eigenvalues.  Then 
X(A11  - A12D  5*  A (A^  - A^l*’)  implies  L ^ L' , so  these 
solutions  are  distinct. 


CHAPTER  4 


COMPUTING  THE  LK  TRANSFORMATION 


If  efficient  eigenanalysis  programs  are  available,  a 
computational  procedure  for  obtaining  accurate  numerical  ap- 
proximations to  the  L and  K matrices  is  quite  straight- 
forward. One  such  program,  EISPACK,  is  available  from  the 
Argonne  National  Laboratory  and  is  described  by  Smith  et  al. 
(1976) . This  program  consists  of  many  modular  subroutines 
which  can  handle  a variety  of  eigenanalysis  tasks.  It  is 
well  documented  and  easily  implemented  on  either  batch  or 
timesharing  computers.  Typical  execution  times  for  finding 
all  eigenvalues  and  eigenvectors  of  matrices  of  order  40 
and  80  are  .66  and  4.6  seconds  for  an  IBM  370/195 
computer . 

4.1.  Computing  the  L Matrix 
Since  the  L matrix  can  be  expressed  as  either 


LM 


11 


= -M 


21 


(1.1) 


or 


Q22L  = °21 


(1.2) 
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i 


65 


not  all  eigenvectors  of  A need  be  obtained  to  solve  for  L. 
Two  procedures  could  be  used,  one  for  the  case  n^  <_  n 2 and 
another  for  n-j^  > n2>  For  n.^  <_  n2,  obtain  n1  (right) 
eigenvectors 


M 

M 


11 

21 


corresponding  to  the  slow  eigenvalues  S and  solve  the  set 

of  linear  equations  (1.1)  for  the  L matrix.  Such  equations 

can  be  efficiently  solved  by  a Gaussian  elimination  scheme 

such  as  the  LU  decomposition,  cf.  Stewart  (1973).  In  the 

event  that  S contains  one  or  more  complex  eigenvalue  pairs 

A.,  A.  with  corresponding  complex  eigenvectors  v , v-r- 
ii  Ai  Ai 

replace  v.  and  Vr-  with  Re  (v.  ) and  Im(v.  ).  Then 
Ai  Ai  Ai  Ai 

assuming  the  states  are  ordered  so  that  the  resulting 
is  of  full  rank,  solve  (1.1).  A method  of  approximately 
minimizing  the  condition  number  of  M.^,  avoiding  a singular 
matrix  and  obtaining  an  L matrix  of  small  norm  is 

presented  in  Section  4.3. 

For  n.^  > n2,  an  alternate  procedure  could  be  fol- 
lowed yielding  fewer  computations.  That  is,  find  the  n2 
left  eigenvectors  [Q21  of  A anc*'  after  eliminating 

imaginary  parts,  solve  (1.2)  for  L.  Here  the  term  left 
eigenvector  is  well  defined  by  the  notation  and  equals  the 


66 


corresponding  eigenvector  of  A . This  terminology  differs 
somewhat  from  that  of  some  texts  which  refer  to  the  rows  of 
Q as  the  left  eigenvectors  of  A. 

Let  Lq  be  the  matrix  obtained  by  approximately 
solving  either  (1.1)  or  (1.2).  The  accuracy  of  Lq  can  be 
checked  by  substituting  it  back  into  the  nonsymmetric  alge- 
braic Riccati  equation  and  evaluating  the  residual  error 
matrix 


R 

o 


L Ani 
o 11 


A22Lo 


L A..L 
o 12  o 


+ A 


21' 


J 


(1.3) 


If  Rq  is  judged  to  be  not  sufficiently  small,  i.e. , if 
|| R0 II  > < for  some  small  tolerance  k > 0 and  some  comput- 
able norm  ||*||,  then  the  accuracy  of  Lq  can  be  improved  by 
an  iterative  procedure. 

An  efficient  method  of  improving  the  initial  approxi- 
mation Lq  to  the  decoupling  matrix  L can  be  obtained  by 
using  a successive  approximations  scheme  on  ARE  with  Lq  as 
the  first  iterate.  Using  the  fact  that  X (A^2)  = F we  can 
write  ARE  as 

L = * A21>  = <A22  + LA12)~1(LA11  + A21>- 

If  Lq  is  a reasonable  approximation  to  L,  then  the  matrix 
(A^2  + L0Ai2^  ^as  lar9e  eigenvalues,  is  nonsingular,  and 


I 


the  preceding  equation  suggests  the  linear  iteration 


Li+1  (A22  + LiA22)  (LiAll  + A21J  * 


This  iterative  correction  takes  specific  advantage  of  the 
two-time-scale  property  of  the  system  and  can  be  implemented 
as 

Algorithm  4.1. 

1)  Obtain  the  initial  approximation  Lq  from  either  (1.1) 
or  (1.2)  and  set  i = 0. 

2)  Evaluate 


Ri  LiAll  A22Li  ” LiA12Li  + A21 


and  stop  if  ||R^||  < k. 

3)  Solve  (A22  + LiA12)Di  = ^ for  Di  and  let 

L . . = L . + D . • 
l+l  1 l 

4)  Let  i = i + 1 and  go  to  2. 

It  can  be  shown  that  for  close  to  L,  the  proposed 

iteration  converges  with  approximate  rate  of  convergence  p , 
This  follows  since 


Ri+1  * <4  + Di)A11  - A22(Ll  * Dl) 


' tti  + Di>A12<Ll  + Di>  + A21 


68 


or 


Ri+1  " Ri  + Di(All  ~ A12Li)  " (A22  + LiA12)Di  “ DiA12Di' 


Taking  the  correction  matrix  as  in  step  3, 


Ri+1  ~ (A22  + LiA12’  Ri(All  " A12Li  ” A12Di' 


" (A22  + 1iA12,'lRi(All  - A12Li+l’-  (1'4) 


The  largest  eigenvalues  in  absolute  value  of  (A22  + LA12^ 


-1 


and  (An  - A L)  are  and  s , so  that  the  approxi- 

xx  x x ^ 

mate  rate  of  decrease  in  norm  of  the  residual  error  matrix 
in  Algorithm  4.1  is 


nn 


lRi« 


- If 


= v <<  1. 


(1.5) 


This  approximate  rate  of  convergence  can  be  further 
obtained  as  follows.  Let  ||  • ||  be  some  consistent  matrix 
norm,  and  let  differ  from  L by  order  e for  some 

0 < e <<  1,  i.e.. 


E±  = L^  - L = 0(e). 


Then  equation  (1.4)  can  be  written  as 
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Ri+1  (A22  + EiA12}  Ri(All  ' A12Ei) 


Thus 


Ri+1  = {^22~1  + + 0(e)) 


- X22'lRAl  * 0U)- 


Expressing  A.^  and  &22  * in  spectral  form 


% -1^-1  -1  -1 
A”  “ * A22  = Q22  J2  Q22 


11  llulll  ' 


where  J..  = diag(s1,  . ..,  s ),  J = diag(f..,  . ..,  f 
^ -i-  z i n 

and  pre-  and  post-multiplying  (1.7)  by  Q22  and  M^, 

obtain 


Q22Ri+lMll  * J2'lQ22RiMllJl  + 0U)- 


If  the  A matrix  is  nondefective  then 


llJ2“1H  " O and  llJilt  = |si  | 


(1.6) 


(1.7) 


) , 

we 


(1.8) 


for  any  of  the  standard  matrix  norms,  so  taking  norms  in 
(1.8)  it  follows  that 


I1Q22W11I1 

^Q22RiMll^ 


< V + 


0(e)  . 
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(1.9) 


Since  Q22  and  M^l  are  constant  and  nonsingular,  (1.9) 
implies  (1.5)  by  the  equivalence  of  matrix  norms.  Also  the 
residual  matrix  and  error  matrix  can  be  related  as 


R. 

1 


LiAll 


- a22l. 


L . A, „L. 
l 12  1 


+ A 


21 


EiAll  " ^22Ei  " EiA12Ei 


so  that  for  e sufficiently  small  a decrease  in  R^  yields 

a proportional  decrease  in  . 

The  convergence  of  Algorithm  4.1  might  be  improved 

by  including  a line  search  in  step  3.  That  is,  if  ||R^+1|| 

> || R . ||  then  let  L.  , = L.  + hD.  and  search  on  h until 
i"  l+l  1 1 

|| Rj__j_j_ II  minimized.  However,  this  safeguard  procedure  was 

not  needed  in  the  numerical  examples  which  follow. 

Our  procedure  for  obtaining  numerical  solutions  to 
the  nonsymmetric  Riccati  equation  has  some  similarities  to 
one  proposed  by  Farrar  and' DePietro  (1977)  for  solving  the 
symmetric  Riccati  equation  in  that  both  procedures  are 
initialized  using  eigenvectors.  However,  their  iterative 
correction  involves  a bilinear  equation  whereas  Algorithm 
4.1  involves  a linear  iterative  correction. 
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4.2.  Computing  the  K Matrix 
Since  the  Lyapunov  equation 


kX22  - alxK  + a12  = 0 


(2.1) 


has  a unique  solution,  the  choice  of  solution  method  is  not 
so  important.  The  two- time- scale  structure  may  again  be  used 
to  good  advantage,  as  it  was  in  solving  for  the  L matrix  in 
the  following 


Algorithm  4.2. 

1)  Set  Kq  =0,  = ~A^2>  and  i = 0. 

2)  set  Ki+1  = K.  + ftjija-1. 

3)  Evaluate  " Ki+i^22  ~ A12  ant^  stoP  if 

ll^i+lll  < k. 

4)  Set  i = i + 1 and  go  to  step  2. 

An  advantage  of  this  approach  is  that  ft^2  ^ is  un_ 
changed  in  the  iteration  so  that  it  can  be  decomposed  once 
into  LU  form  and  then  used  repeatedly  in  step  2 to  find 

- Kiti  - V 

The  approximate  rate  of  convergence  can  then  be 
demonstrated  as  it  was  for  Algorithm  4.1.  Thus, 


fti+l  - lKi  + 6i>*22  - Xll<Ki  + 6i>  + *12 


= ft  + ft  ft  -ft  fa  = _ft  ft  ft  -1 
Ri  1 22  AllDi  AllRiA22 
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so 


S1 

llRill 

f 

nl 

Hence  the  algorithm  is  simple  to  implement  and  converges 
rapidly  for  systems  that  are  strongly  two- time-scale . 

4.3.  Ordering  and  Partitioning 
of  the  State  Variables 

In  Chapter  2,  the  vector  of  state  variables,  x,  was 
partitioned  into  sub-vectors  and  x ^ with  dimensions 

determined  by  the  number  of  slow  modes,  n^.  However,  there 
may  be  more  than  one  possible  choice  for  n^.  In  some  cases, 
all  eigenvalues  may  be  well  separated  in  absolute  values  so 
any  1 £ n^  < n will  yield  a reasonably  small  parameter  y. 

In  most  cases,  the  choice  of  n^  will  be  strongly 

guided  by  practical  considerations.  If  one  is  working  with 

systems  of  large  order  and  is  seeking  a reduced-order  model 

of  the  slow  mode  dynamics,  the  choice  of  n^  would  be  based 

on  a trade-off  between  model  accuracy  and  computing  time  and 

expense.  As  described  earlier,  the  number  of  multiplications 

per  time  step  for  the  reduced  model  increases  at  least  as 
2 

rapidly  as  n^.  Also,  shorter  time  steps  are  required  to 
integrate  faster  dynamic  modes.  The  advantages  of  keeping 
n^  small  are  thus  quite  apparent. 
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r 

) 


The  choice  of  might  be  established  by  the  need 

to  model  a given  dynamic  bandwidth  for  control  system  design 
studies,  (e.g.,  Skira  and  DeHoff,  1977).  If  one  needs  to 
design  a controller  to  relocate  n'  slow  open-loop  eigen- 
values, then  provided  the  two-time-scale  assumption  is 
satisfied,  any  n^  n'  can  be  chosen. 

After  the  number  n^  has  been  chosen  the  question  of 
which  state  variables  x should  be  placed  in  arises, 

i.e.,  should  x be  reordered  before  it  is  partitioned? 

As  previously  shown,  the  only  fixed  requirement  on 
the  ordering  of  x is  that  the  resulting  submatrix  of 

the  modal  matrix  M be  nonsingular.  If  the  state  variables 
are  reordered,  then  the  rows  of  the  modal  matrix  M are  re- 
ordered in  the  same  manner.  If  the  purpose  of  decoupling  is 
to  obtain  a reduced  order  model,  this  provides  no  firm  con- 
straints on  the  ordering  of  x since  the  model  yields  the 
slow. mode  dynamics  :ofr.all -state  variables.  There  are,  how- 
ever, sound  reasons  for  considering  a reordering  of  x be- 
fore partitioning.  .If. matrices  L and  K can  be  made  small 
in  norm,  then  both  the  transformation  from  state  coordinates 
x to  time  scale  ecoupled  coordinates  y defined  by  y = Tx 
and  the  inverse  transformation  back  to  state  coordinates 
x = ^y  can  be  performed  with  minimal  numerical  error. 


4.4.  The  Relationship  Between 
Ordering  and  Coupling 
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The  relative  size  of  ||M^^|[  and  llM2 1 H Provi^es 
some  measure  of  the  distribution  of  the  slow  modes  between 
state  variables  x^  and  x From  our  result  (1.1)  it  is 
apparent  that  ||L||  is  proportional  to  llM2 1 1*  an<^  inversely 
proportional  to  If  one  chooses  an  ordering  of  x 

which  minimizes  relative  to  |(M^jJ|  this  will  insure 

that  the  slow  modes  appear  primarily  in  the  components  of  x^, 
that  the  linear  equation  (1.1)  is  well-conditioned,  and  ||L|j 
will  be  small.  In  a similar  fashion,  since 


K = -M  O = -M  (M  - M M ^ 

12y22  12u  22  21  11  12;  ' 


(4.1) 


choosing  an  ordering  of  states  which  minimizes  II 2 ^ w-*-H 

approximately  place  fast  variables  in  and  yield  small 

11*11- 

This  process  is  clarified  by  introducing  slow  and 
fast  mode  coupling  ratios  pg  and  p^.  Assume  that  the 
columns  of  M are  normalized  so  that  each  has  Euclidean 
norm  one,  and  define 

ll«21ll  ii«12ii 

°s H p^r  p£ = p^it- 

• mt*  i • f is  a much  simpler  problem  than  minimizing 

■ n<i  has  the  advantage  of  specifically 
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placing  states  which  are  predominately  fast  into  x2 , cf . , 

(II. 3. 6).  Either  pg  or  can  be  minimized  by  searching 

over  all  orderings  of  the  row  indices  {i}  of  the  modal 

matrix  M = [m. . ] . 

13 

If  x can  be  ordered  so  that  = 0,  then  pg  = 0 

and  the  slow  modes  are  naturally  decoupled  from  the  x2 
states.  For  this  ordering  of  x the  A matrix  will  be 
block- triangular . Likewise,  if  M^2  = 0,  then  p^  is  zero 
and  the  resulting  A will  be  block- triangular . Depending 
upon  the  structure  of  M,  the  ordering  of  x which  mini- 
mizes pg  may  or  may  not  also  minize  p^. 

If  both  p and  p£  are  small  we  could  describe 
s f 

the  system  as  being  "weakly  coupled" . This  description  of 

weakly  coupled  systems  is  similar  to  one  proposed  by  Milne 

(1965) . However,  his  definition  also  required  that  p be 

small.  Even  though  a system  is  strongly  two-time-scale  it 

need  not  be  weakly  coupled,  or  conversely. 

The  following  procedure  is  suggested  for  reindexing 

x.  A,  B,  C,  and  M in  order  to  approximately  minimze  ||L[|, 

place  slow  variables  in  x^,  and  avoid  a singular  M^. 

Similar  procedures  could  be  adopted  in  order  to  approximately 

minimize  either  p£  or  max{p  , p-}.  It  is  assumed  one  is 

f s f 

able  to  compute  approximate  eigenvectors. 

Algorithm  4.3. 

1)  Normalize  the  columns  of 


I 


M 


11 


M 


21 


76 


(4.10) 


n 


so  that  l |m  . | =1,  j = 1,  n.  . 

i=l  13  1 

2)  If  (4.10)  contains  one  or  more  complex  eigenvectors  v, 
replace  v with  Re(v)  and  replace  v with  Im(v)  (as 
described  earlier) . 

3)  Evaluate  the  norm  of  each  of  the  n rows  of  the  re- 
sulting matrix  (4.10)  and  let 


r . = 

i 


n 1 
V 2 
.1 

j-1  1 


1/2 


1 X/  • • • f n a 


and 


4)  Order  the  {r . > so  that  r > r >•...>  r 

i a.  — a_  — — a ' 

12  n 

premultiply  x.  A,  B,  and  M by 

fe~  . . . e — ! T 

L “l  “2  an  I 


where  is  the  ith  column  of  the  n x n identity 

matrix. 

5)  Post-multiply  A and  C by  the  inverse  matrix 


[%  S - e%]‘ 
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One  advantage  to  such  a systematic  approach  to  ap- 
proximately minimizing  ||L||  and/or  ||K||  is  that  it  can  be 
readily  implemented  for  systems  of  large  order,  where  the 
modal  matrix  becomes  unwieldy  for  manual  manipulations. 


4.5.  The  Computational  Method  Summarized 
To  summarize,  the  proposed  procedure  for  computing 
the  decoupling  transformation  for  a two-time-scale  system 
involves  the  following  steps: 

1)  Obtain  approximate  eigenvalues  A (A)  and  choose  the 
number  of  slow  modes  n^.  Note  that  different  choices 
of  n^  are  possible  to  suit  the  needs  of  the  study. 

2)  If  n^  £ n2  obtain  approximate  eigenvectors 


M 
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Otherwise,  obtain  left  eigenvectors  [Q21  Q22]  . 


3)  (Optional.)  Apply  Algorithm  4.3  to  avoid  singular 
and  place  slow  variables  in  x^. 

4)  Eliminate  imaginary  parts  as  needed  and  approximately 
solve  either  LM 


11 


-»21  or  Q22L  " °21 


for  L = L . 

o 


5)  (Optional.)  Improve  the  accuracy  of  Lq  by  Algorithm 
4.1. 

6)  Find  K by  Algorithm  4.2. 


Numerical  examples  of  this  procedure  are  presented  in 


the  next  chapter. 
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NUMERICAL  EXAMPLES 

Two  examples  of  linear  control  systems  taken  from 
recent  literature  are  presented  here  to  demonstrate  actual 
computation  times  and  convergence  rates  of  the  proposed  algo- 
rithms and  also  provide  further  comparison  of  the  new  reduced 
order  model  with  the  Marshall  model. 

The  numerical  computations  were  performed  at  the 
University  of  Arizona  Computing  Center  on  both  the  DEC10 
and  CYBER175  computers  to  compare  the  effects  of  computer 
word  length  and  compare  execution  times  for  time  sharing 
against  batch  processing. 

The  performance  of  Algorithms  4.1  and  4.2  is  measured 
by  the  parameter  log^g  II II  where  R^  is  the  appropriate 
residual  error  matrix  and  ||*||  is  the  Euclidean  norm.  As 
illustrated  below,  the  rate  of  convergence  of  these  algo- 
rithms is  well  approximated  by  the  system's  small  parameter, 
y.  The  numerical  data  for  these  examples  is  listed  in 
Appendix  A. 

5.1.  A Strongly  Two-Time-Scale  Example 

The  first  example  is  a fourth  order  model  of  the 
longitudinal  dynamics  of  an  F-8  aircraft  used  by  Teneketzis 
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and  Sandell  (1977)  to  study  linear  regulator  design  by  a 
multiple-time-scales  method.  Their  design  method  is  based 
on  the  singularly  perturbed  form  of  the  system  described  by 
(II. 3. 16),  and  the  size  of  the  small  parameter  was  not 
evaluated.  The  F-8  aircraft  was  used  by  NASA  as  an  experi- 
mental vehicle  for  a digital  fly-by-wire  technology  research 
program  conducted  from  1973  to  1978  as  described  by  Elliot 
(1977) . The  four  state  variables  for  this  model  are  velocity 
variation  v,  flight  path  angle  y,  angle  of  attack  a, 
and  pitch  rate  q;  and  the  one  control  variable  is  elevator 
angle  u.  The  system  eigenvalues  are  s^,  s2  = “*0075 
+ i.076  and  f^,  f 2 = -.94  + i2.98  so  the  time  scales 
ratio  | S2 | / IfjJ  equals  .024.  Since  this  value  of  y is 
very  small  relative  to  one,  the  system  is  strongly  two-time- 
scale, and  the  large  difference  between  speeds  of  the  slow 
and  fast  modes  is  clearly  seen  in  the  response  to  initial 
conditions  as  shown  in  the  illustrations  which  follow. 

This  system  provides  an  example  of  highly  oscillatory 
dynamics  since  the  damping  ratios  for  the  slow  and  fast  modes 

are  only  .1  and  .3,  respectively.  The  minimum  slow  mode 
. . -4 

coupling  ratio  p is  2 x 10  indicating  that  the  slow 

s 

mode  appears  predominately  in  two  variables,  namely  velocity 
variation  and  flight  path  angle.  These  slow  variables  are 
the  first  two  states  in  vector  x,  so  state  reordering 
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described  in  Chapter  4 was  not  needed,  i.e..  Algorithm  4.3 
was  not  employed. 

The  convergence  of  Algorithms  4.1  and  4.2  for  this 
example  is  illustrated  in  Figure  5.1.  As  shown,  the  initial 
matrix  Lq  obtained  using  EISPACK  is  quite  accurate  and  only 
two  iterations  of  Algorithm  4.1  are  needed  to  reach  the  limit 
accuracy  imposed  by  the  finite  word  length  of  the  computer. 
The  six  orders-of-magnitude  difference  between  the  minimum 
residual  errors  obtained  with  the  DEC10  and  the  CYBER  can 
be  directly  attributed  to  their  36  bit  and  64  bit  word 
lengths.  Computation  times  for  this  example  were  too  short 
to  provide  a meaningful  comparison,  but  computation  times 
are  later  tabulated  for  the  second  example. 

The  new  reduced  order  model  and  Marshall's  reduced 
order  model  for  the  example  are  compared  in  Figures  5.2,  5.3 
and  5.4.  The  response  to  initial  conditions  for  a slow  vari- 
able, velocity  variation,  is  illustrated  in  Figure  5.2  and 
two  time-scales  are  used  to  show  the  slow  and  fast  dynamic 
parts.  As  indicated  in  the  figure,  the  new  reduced  order 
model  provides  a more  accurate  representation  of  the  vehicle 
response  on  both  the  slow  and  fast  time  scales  except  for 
only  the  first  .05  sec.,  i.e.,  the  true  response,  equals  the 
value  for  Marshall's  model  at  time  zero  and  then  rapidly  de- 
parts from  that  value  and  approaches  the  new  model  values. 

The  response  to  initial  conditions  for  a fast  variable,  angle 


angle  of 
(rad) 


angle  of 
(rad) 


Figure  5. 


attack 


attack 


3.  Response  of  a fast  variable,  angle  of  attack, 
of  Example  5.1  to  initial  conditions  v(0)  = 
ft/sec,  q(0)  = .05  rad/sec,  a(0)  = 0, 
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angle  of 
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variation 


II  = Marshall's  model 
III  = New  model 


attack 


Figure  5.4.  Response  of  Example  5.1  to  a .1  (rad)  elevator 
deflection  applied  for  one  second. 
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of  attack,  is  compared  in  Figure  5.3.  For  this  example  both 
of  the  reduced  order  models  accurately  reproduce  the  slow 
\ component  of  the  fast  variable.  Since  Marshall's  model  and 

; the  new  model  integrate  differential  equations  with  the  same 

eigenvalues,  their  responses  to  initial  conditions  have 
identical  period  but  differ  in  phase  and  amplitude. 

The  response  to  control  inputs  of  the  two  reduced 
order  models  and  the  original  model  are  compared  in  Figure 
5.4.  The  fact  that  the  new  reduced  order  model  correctly 
decomposes  control  inputs  into  slow  eigenspace  and  fast 
eigenspace  parts  accounts  for  the  very  good  match  in  velocity 
variation  between  the  true  response  and  the  new  reduced  order 
model. 

5.2.  A Weakly  Two-Time-Scale  Example 
This  example  of  a linear  control  system  was  selected 
as  the  theme  problem  for  the  International  Forum  on  Alterna- 
tives for  Multivariable  Control  held  on  October  13  and  14, 
1977  in  Chicago,  cf . Sain  (1977) . The  linear  model  of  an 
F100  turbofan  engine  was  chosen  as  a realistic  example  of  a 
modern  control  problem  involving  several  control  variables 
which  are  varied  simultaneously  and  in  a coordinated  manner 
to  achieve,  in  this  case,  the  goal  of  rapid  thrust  response 
without  violating  any  of  several  operating  constraints.  The 
engine  model  is  sixteenth  order  and  was  obtained  from  a 
detailed  nonlinear  digital  simulation  of  the  engine  using  a 
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numerical  linearization  technique  described  by  Hackney, 

Miller  and  Small  (1977) . This  linear  model  represents  the 
turbofan  operating  at  sea  level  with  a power  level  angle  of 
83  degrees,  i.e.,  rear  maximum  non-afterburner  power. 

Actually  the  linear  model  represents  only  one  of  36  different 
linear  models  corresponding  to  various  altitude,  Mach  number, 
and  power  settings  which  represent  conditions  within  the 
flight  envelope.  The  operating  point  for  this  linear  model 
represents  a condition  which  every  engine  must  pass  through 
on  takeoff. 

The  sixteen  state  variables  for  the  turbofan  consist 
of  two  shaft  speeds,  three  pressures,  and  eleven  temperatures. 
Five  of  the  state  variables  can  be  sensed  including  two  shaft 
speeds,  two  pressures,  and  one  temperature.  The  model  in- 
cludes five  control  variables  of  fuel  flow,  exit  nozzle  area, 
two  variable  vane  angles,  and  compressor  bleed;  and  there  are 
five  output  variables  including  net  thrust,  total  engine  air- 
flow, temperature,  and  two  stall  margins.  This  model  is  not 
used  here  as  a control  design  problem,  but  rather  to  provide 
a realistic  evaluation  of  Algorithms  4.1  and  4.2. 

The  eigenvalues  for  this  example  are  -.65,  -1.90, 

-2.62,  -6.72  + il. 31,  -17.8  + i4.8,  -18.6,  -21.3  + i0.8, 
-38.7,  -47.1,  -50.7,  -59.2,  -175.7,  and  -577.0.  If  the 
number  of  slow  modes  is  chosen  as  15,  5 and  3,  the 
resulting  small  parameters  are  0.304,  0.371,  and  0.383. 


Since  these  values  are  not  particularly  small  relative  to 
one,  we  call  this  system  weakly  two-time-scale. 

The  performance  of  Algorithms  4.1  and  4.2  for  the  two 
cases  n^  = 3 and  n^  = 5 with  the  state  variables  in  their 
original  order  are  illustrated  in  Figures  5.5  and  5.6.  As 
shown,  the  convergence  rate  approximately  equals  the  small 
parameter.  Again,  the  initial  matrix  Lq  obtained  with 
EISPACK  is  highly  accurate;  the  limit  of  numerical  accuracy 
is  reached  in  6-8  iterations  of  Algorithm  4.1,  and  the 
effects  of  computer  word  length  are  readily  apparent.  The 
computation  times  corresponding  to  Figures  5.5  and  5.6  are 
listed  in  Table  5.1.  As  expected,  the  CYBER175  required 
considerably  less  central  processor  time  to  perform  these 
computations  than  the  DEC10. 

The  effects  of  reordering  the  state  variables  before 

T T 

partitioning  x into  (x^,X2)  are  illustrated  in  Figure  5.7 
for  the  case  n^  = 5.  As  shown,  the.  choice  of  ordering  does 
not  significantly  affect  the  rate  of  convergence  of  Algorithm 
4.1.  However,  the  norms  of  the  L and  K matrices  and  the 
values  of  pg  and  pf  are  influenced  by  the  choice  of 
ordering  as  shown  by  Table  5.2.  Here  less  accurate  eigen- 
vectors were  used  to  calculate  L so  the  rate  of  con- 

o 

vergence  of  Algorithm  4.1  could  be  more  clearly  observed. 

The  computation  of  the  L matrix  using  alternative 

initial  values  for  L is  described  in  Figure  5.7. 
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Algorithm  4.1  iteration,  i 


Algorithm  4.2  iteration,  i 


Figure  5.5.  Convergence  of  Algorithm  4.1  and  4.2  for 
Example  5.2  with  n1  = 3,  y = .383. 
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Figure  5.6. 


Algorithm  4.2  iteration,  i 

Convergence  of  Algorithms  4 . 1 and  4 . 2 for 
Example  5.2  with  n.  = 5,  u = .371 
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Table  5.1.  Execution  Times  Corresponding  to  Figures  5.5 
and  5.6. 


Time  to  compute  Lq 

i 

.109 

.615 

.095 

.599 

Time  per  Li  iteration 

.023 

.187 

.021 

.171 

(Algorithm  4.1) 

Time  per  iteration 

.015 

.126 

.012 

.091 

(Algorithm  4.2) 

Computer 
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Table  5.2.  Data  for  Cases  I,  II  and  III. 


Case 

indices 

ps 

pf 

l|K|| 

I 

.260 

342. 

.536 

3.42 

II 

1,2,5,6,10 

.248 

213. 

.340 

1.55 

III 

1,2,9,11,14 

.183 

39.1 

.693 

17.1 

Case  III  was  obtained  by  Algorithm  4.3 
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Algorithm  4.1  is  quite  stable  for  this  example  and  succeeds 
in  converging  to  the  correct  solution  even  when  initialized 
with  Lq  = 0 or  Lq  = A22  ^A21‘  t^^-s  last  approximation  to 

the  L matrix  was  proposed  by  Kokotovic  (1975),  cf. 

(III. 2. 3a)  However,  the  Lq  matrix  obtained  using  eigen- 
vectors provides  a much  more  reliable  initial  value  for 
Algorithm  4.1.  Nonetheless,  good  convergence  can  be  expected 
from  eigenvectors  which  are  not  highly  accurate. 

An  effort  was  made  to  compare  the  new  reduced  order 
model  with  Marshall's  model  for  this  example.  However  due 
to  the  high  order  and  stiffness  of  this  example  a reference 
time  history  was  not  obtainable  using  standard  integration 
methods . 


CHAPTER  6 


EXTENSIONS  TO  TIME-VARYING 
CONTROL  SYSTEMS 

The  time-scale  decoupling  methods  presented  here  may 
be  applied  to  control  systems  where  all  matrices  are  functions 
of  time.  We  now  assume  the  control  system  has  the  form 

x = A ( t)  x + B (t)  u (0.1) 

with  A(t),  B ( t)  and  u(t)  being  specified  continuous  func- 
tions of  time  on  [0,TJ  and  with  the  initial  state  x(0) 
specified. 


6.1.  The  Time-Varying 
LK  Transformation 

The  control  system  (0.1)  will  be  classified  as  two- 
time-scale  on  the  interval  [0,T]  if  there  exist  integers 
n^  and  n2  = n - n-^  such  that  the  time-varying  spectrum 
X (A  (t) ) can  be  partitioned  into  sets  of  n.^  and  n2  members 

X (A  ( t)  ) = S(t)  U F ( t)  on  [0  ,T]  , 


such  that 
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max  | s . (t)  | = £(t)  <<  v.(t)  = min  jf  . (t)  I (1.1) 

s.€S  1 f.€F  3 

1 3 


throughout  0 < t < T and 


max 

0<t<T 


Wi 

v (t) 


<< 


1. 


It  is  not  yet  assumed  that  (0.1)  is  two-time-scale . 

O'Malley  (1969)  and  Harris  (1973)  have  described  a 
time-varying  version  of  the  LK  transformation  for  homoge- 
neous singularly  perturbed  boundary  value  problems  of  the 
form 


Z1  = Cll^zl  + ci2  ^ z2 
e z2  C21(t)zl  + C22(t)z2 


(1.3) 


on  the  interval  [0,TJ  where  z^  and  z 2 are  real  vectors 
of  dimensions  m1  and  m2  respectively,  and  ^ (t) , i, 
j = 1,  2,  are  compatible  matrix  functions  of  time.  The  ap- 
proximate boundary  conditions  are 


D-^fO)  + D2z(T)  = c (1.4) 


, and  e <<  1 is  the  small  parameter.  As 
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shown  by  O'Malley  (1969) , if  the  matrices 

C22(t)  and  Ci;L(t)  - C12(t)C22(t)C21(t) 

are  nonsingular  on  [0,T] , then  there  exists  a nonsingular 
transformation  of  variables 


w = P ( t) z 


which  transforms  system  (1.3)  into  block-diagonal  form 


"i  - anwi 


€W2  - 2;22W2 


on  [0,'f]  with  boundary  conditions  (1.4)  appropriately  trans- 
formed. Similarly,  Harris  (1973)  showed  that  if  the  eigen- 
values of  matrix  C22(t)  have  nonzero  real  part  on  [0,T] 
then  the  LK  transformation  described  below  exists  on  [0,T] 
and  transforms  (1.3)  to  block-diagonal  form.  Further,  the 
initial  value  problem  has  the  limiting  solution  w2  = 0 if 
all  the  eigenvalues  of  &22(t)  have  negative  real  parts 
everywhere . 

To  extend  these  results  to  control  system  (0.1) 
partition  (0.1)  using  integers  n1  and  n2  as 
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(1.5) 


where  x1 , x_  and  u are  vectors  of  dimensions  n, , n 

x ^ 12 

and  m respectively,  and  all  matrices  are  compatible.  As 
for  the  time-invariant  case,  the  transformation  proceeds  in 
two  steps.  First  define  new  variables  by  a nonsingular 

transformation  £(t) 


so  the  transformed  system  (1.5)  becomes 


y = [£a£  1 + ££  1]y  + £bu  = A ( t)  y + £(t)B(t)u. 


(1.7) 


The  partitions  of  the  A(t)  matrix  are  then 


1 I 

I 


and 
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A22  + ^12 


where  all  matrices  are  functions  of  time.  In  order  for  (1.7) 

A 

to  be  upper  block  triangular,  we  need  A2^(t)  = 0,  so  the 
L ( t)  matrix  must  satisfy  the  Riccati  differential  equation 

L = A22(t)L  - LAu(t)  + LA]_2  ( t)  L - A21(t)  (1.8) 

A 

on  [0,T] . Note,  however,  that  from  the  definition  of  A 
if  L ( t ) ? 0 at  some  point  t'  € [ 0 , T]  then  A(t')  need 
not  be  similar  to  A(t')  and  the  eigenvalues  of  the  block- 
triangular  system  need  not  necessarily  equal  those  of  the 
original  system. 

For  the  second  step  of  the  transformation  define  new 
variables  y^  by  the  nonsingular  transformation 
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T(t)  = 


I + K ( t)  L (t)  K ( tf 


L ( t) 


(1.10) 


so  the  system  (1.7)  is  transformed  to 


y = [TAT  1 + 7T  1]y  + TBu  = ^(t)y  + &u  (1.11) 


where 


and 


A11  A11  AH  " A12L 


^12  * + KA22  “ A11K  + A12 


hi  ' 0 


A22  A22  ” A22  + ^12* 


'Si 


Now  A will  be  block-diagonal  on  [0,T]  if  the  K(t)  matrix 
satisfies  the  linear  differential  equation 


K = ( t) K - KA22(t)  - A12(t)  (1.12) 

on  [0,T] , and  the  original  system  (1.5)  is  reduced  to  the 
decoupled  subsystems  (1.11),  i.e.. 


* 

I 


1 
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*1  - xllyl 

+ i^u 

(1.13a) 

y2  = ^22y2 

+ V 

(1.13b) 

where  & is  partitioned  compatibly  into 
conditions  appropriate  for  (1.13)  are 


The  initial 


' 


y (0)  = T(o)x(0). 


Provided  L(t)  and  K(t)  satisfy  (1.8)  and  (1.12),  the 
solution  x(t)  to  the  original  system  may  be  obtained  by 
integrating  (1.13)  and  transforming  y(t)  as 


x ( t ) = T (t) _1y (t) 


fl  -K(t) 


y (t) . 


-L(t)  I + L ( t)  K ( t) 


Systems  of  linear  differential  equations  related  by  a non- 
singular linear  transformation  such  as  (1.3)  are  said  to  be 
kinematically  similar  (Coppel,  1978) . 

In  order  for  the  transformation  (1.9)  and  thereby 


system  (1.13)  to  be  completely  specified,  appropriate  boundary 
values  for  matrices  L(t)  and  K(t)  need  to  be  chosen.  At 
this  point  we  will  assume  that  the  decoupled  system  (1.13) 
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is  two- time-scale.  Then  proceeding  as  in  the  time  invariant 
case,  L (0)  can  be  chosen  to  initially  separate  time  scales 

f\j 

between  A^  and  A22  so  that 


X(Au(0)  - A12(0)L(0))  = £a  (0)  = S (0)  (1.14a) 

and 

X(A22(0)  + L(0)A12(0))  = F (0)  = F (0)  (1.14b) 

where  A(J^(0))  = S'(O)  U Ir'(O)  , and 


L (0)  = 0. 


(1.15) 


(Then,  we  guarantee  that  the  eigenvalues  of  A(0)  are  pre- 

<\j  , , 

served  m A(0),  i.e.,  the  initial  time-scale-decoupling  is 

preserved  under  the  transformation.)  Let  A(0)  have  the 
spectral  form 


A (0) 


M 


11  M12 


M 


21 


M, 


22 


Q 


11 


21 


where  J1  and  J2  are  Jordan  matrices  with 


- s(0). 


X(J2)  = F (0) 
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and  M ^ = q.  Then  assuming  is  nonsingular  (rearranging 

rows  of  x,  if  necessary)  L(0)  will  satisfy 

L (0)  = -M21M11"1  (1.16) 

and  (1.14)  and  (1.15)  will  both  be  satisfied.  Provided 

i 

T ( t)  B ( t)  u ( t)  = B ( t)  u ( t) 

is  slowly  varying  on  [0,T],  this  choice  of  initial  condi- 
tions on  L will  yield  separation  of  fast  and  slow  dynamics 
in  (1.13)  at  the  initial  time  t = 0 and,  by  assumption, 
this  separation  of  time  scales  in  (1.13)  persists  throughout 
[0,T] . If  this  fails  to  persist,  our  procedure  will  not  be 
applicable . 

One  would  hope  to  be  able  to  solve  for  L(t)  and 
K ( t)  using  some  standard  differential  equation  integrator 
in  a numerically  stable  manner.  To  assess  the  stability  of 
(1.8)  assume  that  a matrix  L(t)  satisfies  (1.8)  with  initial 


condition  (1.16)  and  let 

some 

other 

solution 

L(t) 

satis 

- 

fying  (1.8)  be  displaced 

from 

L ( t) 

at  time 

t'  € 

(0  ,T) 

by 

some  small  amount.  Define  the 

perturbation  matrix 

E(t) 

to 

be 


E ( t)  = L(t)  - L(t) 


i 

! 


I 


i 
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for  t € [t',T]  where  E(t')  is  0(6)  for  some  small 

number  6 > 0.  Then  the  perturbation  matrix  satisfies 


E = A22(L  - L)  - (L  - L)A  x + LA12L  - LA12L 


^22E  " E^ll  + EA12E 


(1.17) 


where  - A^L  and  &22  A22  + LA12 


Since  it  is 


assumed  that  the  decoupled  system  (1.13)  is  two- time -scale 

2 

and  EA^2E  is  initially  0(6  ),  the  perturbation  equation 
(1.17)  can  be  approximated  by 


E J K22E 


and  the  stability  of  (1.17)  will  depend  upon  the  stability  of 
the  fast  subsystem  (1.13b).  Thus  to  insure  the  stability  of 
the  Riccati  differential  equation  (1.8),  it  is  assumed  that 
the  fast  subsystem  (1.13b)  is  stable  on  [0,T] . Sufficient 
conditions  for  the  stability  of  time-varying  linear  systems 
(e.g.,  Davison  1968)  are  not  as  simple  as  those  for  constant 
linear  systems,  and  our  previous  conditions  (11.3.1)  for  fast 
mode  stability  are  no  longer  generally  adequate. 

For  y sufficiently  small,  however,  the  singular 
perturbation  results  cited  show  that  fast  mode  stability  is 
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i t 

1 

f 

A 


guaranteed  provided  the  large  eigenvalues  of  &22  have  large 
negative  real  parts.  The  numerical  stability  of  the  linear 
differential  equation  (1.12)  may  be  assessed  by  a similar 
argument.  Let  matrices  K(t)  and  K(t)  both  satisfy  (1.12) 
and  define  the  perturbation  matrix 


G ( t)  = K(t)  - K ( t ) . 


Then 

6 = ^i;lg  - g&22  (1.18) 

and  by  the  assumption  that  & is  two- time-scale  (1.18)  can 
be  approximated  as 

6 % -($22  • 

If  the  fast  subsystem  (1.13b)  is  stable,  K(t)  can  be  ex- 
pected to  grow  without  bound  as  (1.12)  is  integrated  forward 
in  time.  This  is  demonstrated  in  the  numerical  example  which 
follows. 

However,  the  linear  differential  equation  (1.12)  will 


be  stable  when  integrated  backward  in  time  from  T to  0 
and  one  can  obtain  K(t)  in  this  fashion  by  choosing  an 
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appropriate  final  time  value  K(T).  One  natural  choice  for 
K (T)  is  the  unique  solution  to  the  linear  equation 


o = ^l1(t)k  - k&22  (T)  - a12(,T) 


obtained  by  the  iteration 

K = (^U(T)K  - A12(T)  )A22_1(T) 


so  that 


ft(T)  = 0. 

Based  on  these  findings,  the  following  procedure  for 
decoupling  an  initially  two-time-scale  linear  control  system 
(1.5)  is  suggested. 

1)  At  the  initial  time  t = 0,  apply  the  methods  of  the 
previous  chapters  to  find  the  unique  matrix  L(0) 
satisfying 

0 = A22(0)L(0)  - L(0)Ai;l(0)  + L(0)A12(0)L(0)  - a21(0) 


X(A11(0)  - A12(0)L(0))  = A (S  (0)  ) . 


with 
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2)  Integrate  the  time-varying  Riccati  differential  equation 

L = A22(t)L  - LAn(t)  + LA12(t)L  - A21(t) 

from  t = 0 to  t = T and  save  the  solution  L(t) . 

(We  assume  that  L(t)  remains  bounded  throughout  the 
interval. ) 

3)  At  the  final  time  t = T,  find  the  unique  matrix  K(T) 
which  satisfies  the  algebraic  equation 

K (t)  = (&n(lf)KCf)  - A12(T))&22_1(T)  . 

4)  Integrate  the  time-varying  linear  differential  equation 

K = ^xl(t)K  - K&22(t)  - A12(t) 

backward  in  time  from  t = T to  t = 0 and  save  the 
' solution  K(t). 

5)  Using  L(t)  and  K(t),  integrate  the  slow  and  fast 

* mode  decoupled  subsystems 


yl  = (A11  " A12L)yl  + t(I  + KL)Bi  + KB2]u 

i2  - <»22  + “l2>y2  + (LBx  + B2)u 


i. 
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with  initial  conditions 

y1(0)  = (I  + K(0)L(0)  )x1(0)  + K(0)x2(0) 
y2(0)  = L(0)x1(0)  + x2  (0)  , 

and  transform  back  to  obtain  values  x(t)  as  needed, 
where 


' I “ 

-K(t) 

X ( t ) = 

yx(t)  + 

-L(t)_ 

I + L ( t)  K (t) 

y2(t) 


As  an  alternative  to  5) , because  the  fast  mode  sub- 
system (1.13b)  is  assumed  to  be  stable  on  [ 0 , 1* 3 and  a re- 

A 

duced  order  approximation  x(t)  to  the  time  solution  x(t) 
is  desired,  the  fast  subsystem  variables  y 2 (t)  could  be 
approximated  by  the  algebraic  relationship 

y2(t)  = -(A22  + LA12)_1(LB1  + B2)u(t). 

A 

Then  the  approximate  state  x(t)  is  given  by 


I 


I 

-K(t) 

x(t)  = 

yx(t)  + 

-L(t)_ 

I + L(t)K(t)_ 

y2  (t) 
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with  again  determined  by  numerical  integration. 

6.2.  A Second-Order  Example 
The  decoupling  procedure  described  above  is  illus- 
trated with  the  following  second-order  example. 


• 

X1 

• 

_X2_ 

-5  + 3h(t)  -6  + h (t) 


u 


(2.1) 


for  t € [0,5]  where 


h (t)  = 


/ sin , 0 < t < 4 


, 4 < t < 5 


with  initial  conditions 


x^O)  = 1,  X2  (0 ) = 1 


and  control  input 


u(t)  = 


^sin2(^)  , 0 < t < 4 


4 < t < 5, 
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The  initial  eigenvalues  for  this  example  are  A-^,  \^  = -1, 

-5.  They  vary  with  time  as  shown  in  Figure  6.1. 

The  Riccati  differential  equation  for  this  example 
is  the  scalar  equation 

L = (-6  + h ( t)  ) L + L - (-5  + 3h  (t)  ) 


with  the  initial  condition 


L(0)  = 1 

providing  L(0)  = 0.  The  L(t)  time  history  from  t = 0 to 
t = 5 is  illustrated  in  Figure  6.2,  and  the  corresponding 
fast  and  slow  system  matrices 

A11  = A11  “ A12L'  A22  = A22  + ^12 

are  shown  in  Figure  6.1  along  with  the  original  system's 

, % % 

eigenvalues  A^  and  As  shown,  A.^  and  A22  closely 


follow  A^  and  A2  here. 


-6  i 


time 


Figure  6.1.  Eigenvalues  of  original  and  block-decoupled 
system  for  Example  6.2. 

i 


Figure  6.2.  L(t)  and  K(t)  time  histories  for 
Example  6.2. 


i 
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The  scalar  differential  equation  for  K(t)  to  be 
integrated  backward  is 

K = -LK  + K (-6  + h + L)  - (-5  + 3h) 

with  final  value 

K (5)  = \ 

corresponding  to  L(5)  « 1.  The  K(t)  trajectory  is  illus- 
trated in  Figure  6-2.  In  order  to  demonstrate  the  insta- 
bility of  the  differential  equation  for  K(t)  when  integrated 
forward  in  time,  a small  portion  of  such  a trajectory  is  also 
shown  in  Figure  6.2. 

Since  the  system  is  linear,  the  superposition  prin- 
ciple applies,  and  the  solution  x(t)  can  be  considered  as 

the  sum  of  two  parts;  the  response  to  initial  conditions 
T 

x(0)  = (1,1)  with  zero  control  input  u(t),  and  the  re- 
sponse to  control  input  u(t)  with  zero  initial  state. 

Figure  6.3  illustrates  the  response  of  the  reconstructed 
states  x^(t),  x2(t)  and  reduced-order  approximations 

A A 

x^(t),  x2(t)  to  initial  conditions,  whereas  Figure  6.4 

illustrates  the  response  of  these  variables  to  the  control 

A 

input.  As  shown,  the  state  x provides  a reasonably  good 
approximation  to  x for  this  example  away  from  an  initial 


state  variables 


< 

r 

[ 
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Figure  6 
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boundary  layer  where  initial  conditions  are  satisfied.  The 
example  was  also  integrated  in  its  original  form  (2.1)  to 
verify  that  the  reconstructed  state  obtained  as  x = T ^ y is 
correct. 

As  demonstrated  by  this  example,  the  time  scale  de- 
coupling transformation  can  be  extended  to  time  varying 
linear  control  systems  in  a straightforward  manner.  The 
method  should  prove  to  be  particularly  useful  for  large  scale 
linear  systems  where 

1)  the  L ( t)  and  K(t)  matrices  are  slowly  time  varying, 
and 

2)  the  decoupled  system  is  strongly  two- time-scale. 

In  this  case  only  the  fast  subsystem  would  require  short  in- 
tegration steps.  Also,  for  strongly  two-time-scale  systems 
where  the  fast  subsystem  is  strongly  stable,  the  reduced- 
order  approximation  proposed  here  could  lead  to  considerable 
reduction  in  computation  time  requirements,  especially  for 
large  order  problems. 


I ' w 


( 

f 
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APPENDIX  A 

DATA  FOR  EXAMPLES  5.1  AND  5.2 

The  state  and  control  variables  for  the  aircraft 
longitudinal  dynamics  example  are 


x = 


u = 


6 

e 


where 

v = horizontal-velocity  deviation  (ft/sec) 

Y = flight-path  angle  (radians) 

a = angle  of  attack  (radians) 

q = pitch  rate  (radians/sec) 

<5  = elevator  deflection  (radians) . 

0 

The  corresponding  system  matrices  are 


-.01357 

-32.2 

-46.3 

0.0 

.0001200 

0.0 

1.214 

0.0 

A = 

-.0001212 

o 

• 

o 

- 1.214 

1.0 

_ .0005700 

0.0 

- 9.01 

- .6696^ 
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and 


-.43301 


B 


.1394 

-.1394 


1_-.1577 


The  sixteen  state  variables  for  the  turbofan  engine 
example  are 
x^  = fan  speed  (rpm) 

*2  ~ compressor  speed  (rpm) 

x^  = compressor  discharge  pressure  (psia) 

x4  = interturbine  volume  pressure  (psia) 

Xj.  = augmentor  pressure  (psia) 

Xg  = fan  inside  diameter  discharge  temperature  (°R) 
x^  = duct  temperature  (°R) 
x8  = compressor  discharge  temperature  (°R) 
xg  = burner  exit  fast  response  temperature  (°R) 
x^g  = burner  exit  slow  response  temperature  (°R) 
x^  = burner  exit  total  temperature  (°R) 
x^2  = fan  turbine  inlet  fast  response  temperature  (°R) 
x13  = fan  turbine  inlet  slow  response  temperature  (°R) 
x14  = fan  turbine  exit  temperature  (°R) 
x15  = duct  exit  temperature  (°R) 
xlfi  = augmentor  temperature  (°R) 


and  the  five  engine  input  variables  are 

u^  = main  burner  fuel  flow  (lb/hr) 

2 

U£  = nozzle  jet  area  (ft  ) 
u^  = inlet  guide  vane  position  (deg) 
u4  = high  variable  stator  position  (deg) 
u5  = customer  compressor  bleed  flow  (%) . 

The  numerical  data  for  the  A and  B matrices  for  the 
turbofan  engine  are  listed  in  Tables  A.l  and  A. 2. 


tsBEsTa'J“'lo5W0!l^M 
^^Uis'®10000  ■ 
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Table  A.l. 


The  A Matrix  for  the  Turbofan  Engine  Model. 


Columns  1 through  8 : 


-4.328 

.1714  5.376 

401.6 

-724.6 

-1.933 

1.020 

-.9820 

-.4402 

-5.643  127.5 

-233.5 

-434.3 

26.59 

2.040 

-2.592 

1.038 

6.073  -165.0 

-4.483 

1049. 

-82.45 

-5.314 

5.097 

.5304 

-.1086  131.3 

-578.3 

102.0 

-9.240 

-1.146 

-2.408 

. 84  76E-02 

- . 1563E-01  . 5602E-01 

1.573 

-10.05 

.1952 

- . 8804F-02 

-.  21 10E-01 

.8350 

- . 1249E-01  -.3567E-01 

-.6074 

37.65 

-19.79 

-.1813 

- . 2962F-01 

.6768 

- . 1 264E-01  - . 9683E-01 

-.3567 

80 . 24 

- . 82398-01 

-20.47 

- . 39288-01 

- . 9696E-01 

.8666  16.87 

1.05] 

-102.3 

29.66 

.5943 

-19.97 

-.8785E-02 

- . 1636E-01  .1847 

.2169 

-8.420 

.7003 

. 5666E-01 

6.623 

-.  1298E-03 

-. 2430E-03  . 2718E-02 

. 3214F.-02 

-.1246 

.10398-01 

.83958-03 

.98128-01 

-1.207 

-6.717  26.26 

12.49 

-1269. 

103.0 

7.480 

36.84 

-. 2730E-01 

-.4539  -52.72 

198.8 

-28.09 

2.243 

.1794 

9.750 

1206E-02 

- . 201 7E-01  -2.343 

8.835 

-1.248 

.9975F.-01 

. 8059F.-02 

.4333 

-.1613 

-.2469  -24.05 

23.38 

146.3 

1.638 

.1  185 

4.486 

- . 1244E-01 

. 3020E-01  -.1198 

-.4821E-01 

5.675 

-.4525 

19.81 

.1249 

-1.653 

1.831  -3.822 

113.4 

341.4 

-27.34 

-2.040 

-.6166 

Columns  9 through  16: 

i 

.9990  1.521 

-4.062 

9.567 

10.08 

-.6017 

-.1112 

.96028-01 

11.32  10.90 

-4.071 

-.57398-01 

-.6063. 

-.74888-01 

-.5936 

-.96028-01 

-. 9389E-02  .1352 

5.638 

. 2246E-01- 

. 1 797  ■ 

.24078-01 

1.10O 

. 2743E-01 

-3.081  -4.529 

5.707 

-.2346 

-2.111 

-.2460 

-.4686 

-.3223 

. 2090E-02  - . 5256E-01 

-.4077E-01 

-.91828-02 

-.81788-01 

.34288-01 

.49958-0? 

-.12568-01 

-. 1953E-01  -.1622 

-.6439E-02 

-.23468-01 

-.220] 

-.25148-01 

-.37498-02 

-.33618-01 

. 1878E-01  -.2129 

-.9337E-02 

-.31448-01 

-.2919 

-.33708-01 

.88738-01 

- . 4 <588-01 

. 2253E-01  .1791 

. 8371E-02 

. 2645E-01 

.2560 

.28358-01 

-.37498-01 

. 3635F-01 

-49.99  . 6760E-01 

39.46 

.49918-02 

.89838-01 

. 53498-02 

.0 

.13728-01 

-.6666  -.6657 

.5847 

.6654E-04 

.13478-02 

.71318-04 

.0 

.20578-03 

.2854  2.332 

-47.65 

.3406 

3.065 

.3624 

-.4343 

.4681 

-9.627  -9.557 

38.48 

-50.01 

.1011 

. 1 2038-01 

-.46868-01 

. 171 5E-01 

-.4278  -.4245 

1.710 

-2.000 

-1.996 

. 53498-03 

- . 1 999F-02 

.75448-03 

-4.414  -4.354 

17.66 

-3.113 

-3.018 

-19.77 

-.49998-01 

.15098-01 

- . 1127E-02  - . 6760E-02 

. 1835E-01 

-.99818-03 

-.13478-01 

-.10708-02 

-20.00 

-.20578-02 

.5004  -.1437 

-2.416 

-.1073 

-1 .078 

30.63 

19.89 

-50.16 

jocoa  ' 

1^01  COPY  FUKW-SH^  iv 
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Table  A. 2.  The  B Matrix  for  the  Turbofan  Engine  Model, 


-.4570F.-01 

-451.6 

-105.8 

-1.506 

851.5 

.1114 

-546.1 

-6.575 

-107.8 

3526. 

.2151 

1362. 

13.46 

20.14 

-.67778+05 

.3262 

208.0 

-2.888 

-1  .653 

-269.1 

. 9948E-02 

-98.39 

. 5069 

-.1940 

-94.70 

.2728E-01 

71.62 

9 . 608 

-.11 60 

-184.1 

. 1716K-01 

71.71 

8.571 

. 798«i 

-515.2 

- . 7741E-01 

-141.2 

-.8215 

39.74 

1 376. 

. 3855E-01 

-7.710 

-.43718-01 

-.1024 

-6684. 

. 5707E-03 

-.1144 

- . 6359E-03 

-.  1432E-02 

-99.02 

5.727 

-1745. 

-8 . 940 

-17.96 

. 8898E+0S 

.1392 

-24.30 

-.2736 

-.3403 

-6931 . 

. 61 72E-02 

-1.082 

-. 1183E-01 

-.  1452F.-01 

-307.7 

. 6777K-01 

16.60 

.3980 

.2311 E-01 

-2588. 

. 1880E-02 

9.147 

-.8241 

. 8984K-01 

-32.31 

.1677 

435.8 

-89.94 

4 . 900 

-295.5 

i 
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